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FOREWORD 


The  work  on  variable  moduli  mo del a  vaa  conducced  for 
Che  U.S.  Army  Engineer  Waterways  Experiment  Station  (WES) 
under  sponsorship  of  the  Defense  Atomic  Support  Agency  (DASA) 
as  part  of  NWER  Subtask  SB209,  "Propagation  of  Ground  Shock 
Through  Earth  Media". 

The  contract  was  monitored  by  Dr.  J.S.  Zelasko,  Impulse 
Loads  Section,  Soil  Dynamics  Branch,  under  the  general  super¬ 
vision  o£  Mr.  J.P.  Sale,  Chief,  Soils  Division,  WES. 

COL  Ernest  D.  Peixotto,  CE,  was  the  contracting  officer. 
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ABSTRACT 

The  present  report  extends  the  combined  veriable  moduli 
model  introduced  in  the  previous  report.  Ref.  [1].  This  more 
general  model  is  defined,  conditions  are  set  on  the  various 
parameters,  and  the  model  behavior  in  uniaxial  strain, 
triaxial  compression  and  proportional  loading  tests  is  dis¬ 
cussed. 

The  major  portion  of  the  present  report  deals  with  the 
procedures  used  to  fit  the  current  model.  Including  the 
loading  and  unloading,  to  a  rather  complete  set  of  laboratory 
data  for  McCormick  Ranch  Sand.  Actually,  four  different  fits 
are  described,  one  of  them,  Uniax-Triax  I,  in  some  detail. 

The  theoretical  and  experimental  results  are  compared  and 
with  one  fit,  Uniax-Triax  II,  excellent  agreement  is  found 
for  uniaxial  strain,  triaxial  compression  and  proportional 
loading  tests. 

Finally,  recommendations  are  made  concerning  reloading 
in  shear.  User's  guides  and  FORTRAN  listings  of  the  two 
programs,  UNAX2  and  PROP,  used  to  compute  uniaxial  strain  and 
proportional  loading  (including  triaxial  compression)  tests 
are  given  in  an  appendix. 
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The  objectives  of  the  work  reported  herein  were  to  (a)  evaluate  the 
capability  of  variable  moduli  type  models,  Ref.  [lj,  to  quantitatively 
match  typical  load-unload-reload  soil  property  test  data  obtained  with 
different  laboratory  test  devices,  and  (b)  document  procedures  for  de¬ 
termining  tiie  model  coefficients  from  a  given  set  of  data.  This  work 
forms  part  of  a  broad  theoretical  and  experimental  research  program  being 
conducted  at  WES  under  DASA  NWER  Subtask  SB209,  "Propagation  of  Ground 
Shock  Through  Earth  Media,"  aimed  at  defining,  describing  and  evaluating 
those  characteristics  of  earth  media  which  govern  the  propagation  and 
attenuation  cf  ground  shock. 

All  laboratory  data  utilized  for  this  report  were  furnished  by  WES. 
Representative  uniaxial  strain  stress-strain  relations  were  determined 
based  on  analyses  of  a  series  of  tests  conducted  at  WES  in  support  of  a 
DASA-sponsored  HEST  test  at  the  McCormick  Ranch  test  site.  Ref.  [2] .  The 
standard  triaxial  compression  test  data,  proportional  loading  shear  data 
and  hydrostatic  compression  data  were  obtained  from  the  Georgia  Institute 
of  Technology,  under  contract  to  WES,  as  part  of  a  research  study  on  the 
behavior  of  soils  under  high  pressure,  Ref.  [3], 

The  vertical  deformation  measurement  of  triaxial  test  specimens  is 
generally  considered  a  routine  operation.  Ref.  [3]  describes  the  develop¬ 
ment  and  pplication  of  an  innovative  device  for  measuring  the  lateral 
deformations  of  cylindrical  soil  specimens  inside  the  triaxial  test 
chamber.  These  additional  measurements  allowed  for  the  computation 


However*  as 


of  the  complete  strain  tensor  for  each  triaxial  test, 
pointed  out  in  Ref.  [3],  and  more  recently  in  Ref.  [4],  interpretation 
of  these  deformation  measurements  i6  still  a  matter  of  major  concern  and 
the  subject  of  intense  additional  research.  Thus,  the  triaxial  test 
data  utilized  in  this  report,  though  generally  self-consistent  and 
assuredly  of  the  correct  order  of  magnitude,  cannot  be  considered 
surgically  precise. 

Finally,  as  pointed  out  by  the  author  in  Section  111  of  this  report, 
there  is  some  disagreement  between  the  virgin  loading  stress-strain 
relations  obtained  with  the  triaxial  test  device  and  the  virgin  loading 
uniaxial  strain  relation.  Considering  the  facts  that  the  two  types  of 
data  were  obtained  with  totally  different  objectives  in  mind,  that  two 
different  pieces  of  test  apparatus  were  used,  each  located  in  a  different 
laboratory  and  each  using  specimens  of  different  size  prepared  by  slightly 
different  techniques,  some  data  disagreement  should  be  no  surprise. 

Indeed,  the  fact  that  the  data  agree  as  well  as  they  do  is  considered 
remarkable  at  WES. 


Vicksburg,  Mies. 
November  1970 


J.  S.  ZELASKO 

Contracting  Officer's  Representative 
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I  INTRODUCTION. 

The  advances  In  Che  design  of  hardened  underground 
structures  have  led  to  Increased  demands  for  accurate  pre¬ 
dictions  of  ground  shock  effects  from  nuclear  explosions. 

These  demands.  In  turn,  have  led  to  the  requirement  of  more 
realistic  mathematical  models  of  the  behavior  of  the  In  situ 
material,  namely,  highly  nonlinear  and  hyateretic  soils. 

The  models.  Ideally,  ahould  reproduce  real  soil  behavior 
for  both  laboratory  tests  and  the  complex  geometries  of  real 
nuclear  and  high  explosive  field  events.  Of  course,  little 
field  data  exists  for  nuclear  events.  The  approach  therefore, 
Is  to  make  the  model  conform  to  material  property  data  obtained 
from  a  variety  of  dynamic  laboratory  tests  and  to  evaluate  the 
model  in  calculations  of  existing  high  explosive  field  tests. 
Although  this  does  not  guarantee  the  correct  behavior  in  a 
nuclear  event,  it.  is  the  most  physically  meaningful  means  of 
obtaining  confidence  In  code  predictions. 

Up  to  now,  all  computations  have  been  performed  with 
elastic-plastic  models.  The  historic  development  of  advanced 
elastic-plastic  models  is  given  In  Ref.  [1].  Although,  as 
Indicated  In  Ref.  [1],  elastic-plastic  models  do  reproduce 
soil  uniaxial  strain  test  data  quite  well,  and  contain  the 
measured  failure  envelope,  they  do  not  correctly  model  material 
behavior  approaching  failure  in  triaxial  compression  tests. 

This  deficiency  led  to  the  development  of  the  variable  moduli 
models  in  which  both  the  bulk  and  shear  moduli  are  functions  of 
the  stress  and/or  strain  invariants,  and  in  which  there  is  no 


—  2  — 

explicit  yield  condition.  The  early  models  discussed  in 
Ref.  [1]  do  give  qualitative  agreement  with  the  results  of 
the  usual  laboratory  tests,  namely,  the  uniaxial  strain  and 
triaxial  compression  tests;  in  this  report,  a  more  advanced 
model  is  described  which  agrees  qualitatively  and  quantitatively 
with  laboratory  test  data  for  a  real  soil,  McCormick  Ranch  Sand. 

In  Section  II,  the  mathematical  c.odel  is  described,  as 
are  analytic  results  for  the  various  test  configurations  and 
allowable  ranges  of  the  material  parameters.  The  available 
laboratory  data  and  the  method  used  to  pick  parameters  for 
loading  and  unloading  are  discussed  in  Section  III.  In 
Section  IV,  results  are  given  for  the  McCormick  Ranch  Sand 
and  the  question  of  reloading  is  reexamined.  Finally,  in 
Section  V,  recommendations  are  given  for  future  work. 
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II  MODEL  DESCRIPTION. 


The  present  Mathematical  model.  Combined  Variable 
Moduli  Model  II,  is  an  extension  of  the  combined  variable 
moduli  model  described  in  Ref.  [1].  As  before,  the  material 
is  described  by  incremental  stress-atrain  relations 


’U 


2Gk 


13 


(1) 


p  "  3K& 


(2) 


The  bulk  modulus  on  virgin  loading  is  retained 

KLD  -  Ko  +  Kle  +  K2e2  (3) 

However,  in  order  to  allow  for  more  general  pressure-volume 
relations  in  unloading,  a  linear  expression  in  pressure  is 
used  as  the  bulk  modulus  for  unloading  and  reloading 


KUN  *  KoU  +  *1UP 

The  major  change,  however,  is  in  the  shear  modulus. 
Whereas  the  failure  envelope  for  the  combined  variable 
moduli  model  described  in  Ref.  [1]  is  a  straight  line 
corresponding  to  a  Prager-Drucker  type  yield  condition, 
the  failure  envelope  of  most  partially  saturated  soils 
starts  as  a  straight  line,  but  then  flattens  out  and 
reaches  a  maximum  with  increasing  pressure.  Advanced 
plastic  models  such  as  the  one  described  in  Ref.  [1] 
mirror  this  behavior  with  a  yield  condition  in  which  fZ 
is  taken  as  a  more  general  function  of  pressure.  Here, 
two  different  expressions  are  used  for  the  shear  modulus. 


For  small  pressure,  i.e., 


for  p  less  than  some  critical 


pressure  p 


G.  -  G  +  Y 
LD  o 


gun  “  GoU  +  Y1U 


(A 


p  +  y2p 


j2  >  0 


*  Y1UP  +  Y2UP  •  j2  i  ° 


13) 


(6) 


If  In  Eq.  (5)  Yj^  >  0  and  Y2  <  °»  then  G^  at  constant  J2 
will  Increase  with  increasing  pressure  until  a  maximum  is 
reached  at 

Y1 

Pc  "  -  W2 

If  it  is  assumed  that  the  same  transition  pressure  p^ 
applies  in  both  loading  and  unloading,  then  necessarily 

llO  IjL 

Y2U  '  Y2 

For  larger  pressures,  i.e.,  p  >  p 


(7) 


(8) 


G,  _  -  G 

LD 


■i +  ?i  pi  • 

fi  > 


gun  "  G1U  +  Y1U 


J2  >  0 


J2  <  0 


(9) 


(10) 


where 


G  -ill 
Go  4  y2 


1  Y1 


1U  oU  4  y. 


(ID 


(12) 


Of  course,  at  p  ■  p{  ,  the  expressions  for  G^  and  G^^  are 
continuous.  It  is  seen  that  the  present  model  has  almost 
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as  nuch  flexibility  in  unloading  as  in  initial  loading, 
something  not  true  of  previous  models.  The  special  case  of 

•  e 

Jj  "  0  la  included  with  Jj  <  0  by  analogy  with  the  neutral 
plastic  state  in  plasticity.  At  present,  the  model  makes  no 
distinction  between  initial  loading  In  shear  and  subsequent 
reloading.  The  ramifications  of  this  and  an  alternative 
choice  are  discussed  later. 

Conditions  on  the  Parameters. 

In  order  for  there  to  be  no  energy  generate.'  during 
infinitesimal  stress  cycles,  two  necessary  conditions  are 


kun  -  kld 

(13) 

gun  1  gld 

(14) 

During  Initial  loading  the  pressure  may  be  found  In  terms 
of  the  mean  strain  by  direct  integration  of  Eqs .  (2)  and  (3) 

p  -  3KQe  +  |  l^e2  +  K2e3  (15) 

Along  the  initial  loading  curve,  by  substituting  Eq.  (15) 
into  Eq.  (4) 

KUN  "  KoU  +  KlU(3Koe  +  \  Kle2  +  K2e3)  (16) 

so  that  the  condition  becomes 

<K„0  -  V  *  <3K10Kc»  *  V  +  <!  KIVK1  -  V*  +  (K10K!)e5  i  0 

(17) 

A  sufficient  condition  for  the  inequality  Eq.  (17)  to  be 
satisfied  for  all  positive  values  of  e  is  for  erch  of  the 
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coefficients  to  be  positive.  The  first  and  second  teras 
lead  to 


K  > 
oU  — 


(18) 

(19) 


If  <  0  and  >  0,  the  third  tern  woul  not  be  satisfied, 
but  since  e  is  snail  this  case  need  not  be  considered. 


The  condition  on  the  shear  moduli,  Eq.  (14),  leads  to 

(GoU  -  V  +  <Y1U  -  *1  <Y10  -  Vp  +  <Y2U  -  Y2)p2  i  0  (20) 

The  initial  shear  modulus  is,  of  course,  positive,  so  that 

GoU  i  Go  >  °  (21) 

Since  in  loading  the  material  softens  with  Increasing  load 
Y^  <  0  and  the  second  term  will  be  positiv  /henever 

Y1U  >  Yx  <  0  (22) 

so  that  Yltj  >  0  is  certainly  possible.  Using  Eqs.  (7)  and 
(8)  the  pressure  terms  in  Eq,  (20)  reduce  to 

p(YlU  "  Yl)(1  "  2?")  -  0  (23) 

v  c 

which  being  valid  only  for  p  <  Pc  leads  to 

Y1U  ”  Y1  ”  °  <24> 

Looking  at  the  loading  and  unloading  shear  moduli  for  p  >  Pc 
leads  again  to  Eqs.  (21)  and  (22)  and  gives  no  new  information. 


7 


Uniaxial  Strain  Test. 

Aa  was  done,  for  the  simpler  combined  variable  moduli 

model  in  Ref.  [1],  the  expression  fer  the  constrained 

4 

modulus  during  initial  loading,  *>  ,  may  be 

Integrated  to  obtain  the  stress-strain  curve,  c^e),  in 
clcsed  form.  The  solution  for  the  stress  is  an  exponential 
plus  a  sixth  degree  polynomial  in  e  with  the  coefficients 
in  terms  of  the  material  parameters.  The  result  is  not 
given  here  (it  would  take  a  whole  page  to  write!)  since  the 
program  which  computed  the  uniaxial  strain  test  actually 
used  numerical  integration  rather  than  the  closed  form  ex¬ 
pressions.  The  other  stress  quantities,  such  as  Sj  and  , 
may  be  obtained  from  C^e)  and  p(e) ,  Eq.  115).  All  strain 
quantities  are  proportional  to  e. 

Trlaxlal  Compression  Test. 

In  the  trlaxlal  compression  test  all  the  stress  rates 
are  proportional  and  since  ■  Oj 


The  shear  modulus  for  p  <  pc  ,  Eq .  (5),  becomes 

—  2 

GLD  ■  Go  +  Ylp  +  ^  Yl(p  "  a3*  +  Y2p  "  Y2(p  ■  pl)(p  ■  p2* 
where  with  Y^^  >  0,  Y-^  and  Y2  <  0»  and  (Yj  +  ^  Yj_)  <  0  as  in 


(25) 


(26) 


I 


—  9 


<01  - 


'r2  ■  0 


3(Co  *  Va> 
H  *  n  *i 


The  expression  for  e^  ,  Eq.  (29) ,  also  reduces  to  that  given 
previously . 


When  o3  <  p,  <  p,  since  the  expression  for  G  changes, 
one  integrates  first  from  <?3  to  pc  and  then  from  p^  to  p. 
The  resulting  expression  for  e^  is 


in  ( 


P„  -  P 


°3  -  pl 


-  In  (-£■ 


P„  -  P- 


°3  -  p2 


) 


y  (Yx  +  /I  Yx) 2  -  ^Y2<g0  -  S3  Y1°3> 


✓3  Yj 


An 


<=.  -  i  r,  *  *  V'  -  °3> 


L  G 


-  i  7: +  71  ^i(pc  ■  °3) 


for  °3  <  Pc  <  P' 


Finally,  when  both  c3  and  p  are  greater  than  p^ 


/3  Yl 


in 


-i^  +  /5  9i(P  - 


°3) 


G  -ill 
Go  4  Y, 


Both  Eqs.  (32)  and  (33)  become  infinite  when 


p  -  p3  -  O,  - 


J  ^1 


(G  -  i  li) 
(Go  4  y. 


(31) 


(32) 


(33) 


(34) 


Em-  (29),  (32)  and  (33) 


The  various  expressions  for  e^  , 
are  continuous  as  either  p  or  0^  crosses  p£  . 


Froa  Eq .  (30)  for  the  stress  difference  at  failure  one 
obtains  tha  intercept 


(01  -  V..I 


'o3  ■  0 


_3 

2y 


(yx  +  /J  rx)  + 

fc-l  *  /It,)2  -  *vT 


and  the  initial  slope 


(35) 


d(0l  “  °3)-.x 


do. 


-3  - 


3/5  y1 


o3  •  o 


/  (Yx  +  /5  Yx)2  -  4Y2Gt 


(36) 


of  the  failure  envelope.  One  may  also  find  the  chaaber 
pressure  o3c  at  which  the  failure  envelope  flattens  out. 


dl"i  -  V...  . 

1 .  e . ,  ■  1  -  •  0 


do. 


3c 


-Y1(Y1  +  2/5  Y^  +  4y2G( 

4y2  yx  /5 


(37) 


It  should  be  noted  that  when  0^  ■  O ^  ,  Eqs.  (28)  and  (34) 
lead  to  p2  ■  P3  ■  Pc  *  It  is  thus  seen  that  in  the  present 
model  there  is  a  unique  pressure  pc  which  represents  the 
pressure  at  which  the  failure  surface  flattens  out,  an 
observed  phenomenon  in  partially  saturated  soils.  The 
present  model  is  therefore  a  strong  candidate  for  representing 
real  soils.  The  actual  fitting  to  laboratory  data  is 
described  in  the  next  section. 
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Proportional  Load  Test?^ 

In  a  so-called  "proportional  load"  test  the  stresses  are 
constrained  according  to  the  relation 


where  q  is  a  constant.  The  two  limiting  conditions  are 
q  ■  0  producing  triaxial  loading 
q  ■  1  producing  hydrostatic  loading 


(38) 


(39) 


The  constitutive  relations  are  those  of  Eqs.  (1)  and  (2), 
the  bulk  moduli  are  given  by  Eqs.  (3)  and  (A),  and  the 
shear  moduli  by  Eqs.  (5),  (6),  (9)  and  (10).  For  pressure 
reloading  up  to  a  previous  peak  K  -  ,  but  for  reloading 

of  deviators  G  •  GLI}.  Figures  (1)  and  (2)  show  typical 
paths  of  loading-unloading-reloading  which  have  been  con¬ 
sidered. 


(1)  Solutions  for  the  Volumetric  Portion. 

Initial  prestressing  of  the  system  has  been  introduced 
through  the  quantities  0^°^  and  °2°^  8uch  that  the  initial 
pressure  is 

p  -  \  (oi0>  +  2 0$o))  (AO) 

ro  3  1  i 

and  initial  mean  strain  eQ  may  be  found  from  the  appropriate 
root  of  the  Eq .  (15)  at  p  -  pQ 

K«e 3  +  1.5K,e2  +  3K  e  -  p  -  0  (Al) 

2o  lo  oo  o 

*7 


This  section  contributed  by  A.  Matthews. 
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The  loading  portion  of  the  volumetric  curve.  Pig.  (1),  from 
* 

Po  to  p  is  fully  described  by  Eq .  (15).  An  unloading 

*  ** 

portion  of  the  volumetric  curve,  from  p  to  p  ,  is  described 
by  the  integral  of  Eq.  (2)  when  K(JN  of  Eq .  (4)  is  used.  The 
expression  is 


e 


3K 


.  K  „  +  K,  ,,p 

1  ln  !_££ - U ZLj 


1U 


oU 


+  Kxup 


* 

+  e 


for  Kly  #  0,  or 


1  /  .  * 
e  “  3 T~  (P  -  p  )  +  e 
oU 

**  * 

for  ■  0.  The  reload  portion  from  p  to  p  is 

** 

established  from  the  same  equations,  except  that  p  and 
**  *  * 

e  replace  p  and  e  ,  respectively,  in  Eqs .  (42)  and  (43). 
For  further  reloading  above  p  ,  Eq.  (15)  is  again  used  to 
describe  the  curve. 


(2)  Solutions  for  the  Devlatoric  Portion. 

Initial  devlatoric  strain,  e^0^  is  considered  to  be 

zero.  With  this  condition,  the  devlatoric  constitutive 

relation,  Eq .  (1),  which  uses  given  by  Eq.  (5)  may  be 

integrated  in  closed  form  for  proportional  loading.  In 

» 

particular,  for  this  case  the  invariant  J2  may  be  written  as 
-V  /  rr  /I  -  q  '  ,  /  _(°)  _(oK 

V  J2  *  (1  +  2q)p  +  1  +  2q  (q0l  "  °2  } 

It  is  therefore  convenient  to  rewrite  the  shear  modulus  G 
in  the  form 

G  -  Gq  +  Gxp  +  G2p2 


(42) 


(43) 


(44) 


(45) 
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where,  for  loading 


G  -  G  +  y  g1  (qa*o) 
o  o  1  +  2q  n  1 

G  m  —  —  ^  /T  y  +  v 
bl  1  +  2q  Vi  Y1  +  T1 


„(oK 
°2  ) 


G  4%  *  Y' 


while,  for  unloading 


5o  -  G0U  +TT^^[0)  -  s20)> 


G1  ■  1+  2q  /3  Y1U  +  Y10 


G2  "  Y2U 


Therefore,  for  p  <  Pc  the  lntegrel  of  combined  Eqs.  (1)  and 


(45)  Is 


•l  -  " 


G1  '  4GoG2 


7  2G2p  + 

.1  2G,p  +  G,  +  *\/i? 


4GoG2 


1  -  AGoG2 


2g2P0  +6,  *  VG1  -  4GqG2 

2G2P0  +  Gx  -V5?  - 


1  o  2 


providing  G2  ^  0.  If  G2  ■  0,  then 


p 

*1  -  ei 


Go  *  Gl* 

Go  +  Vo 


When  q-0  (trlaxlal  case)  Eq.  (48)  la  equivalent  to  Eq.  (29) 


It  may  be  noted  in  Eq.  (48)  that  the  devlatoric  strain 
Increases  without  limit  as  p  approaches  a  value  P^^m  de¬ 
termined  by  setting  the  denominator  of  the  In  term  equal 
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to  zero.  The  resulting  expression  defining  P^ig)  reduces 
to  the  expression  pg  »  Eq.  (28)  when  q-0.  [Alternatively, 
setting  the  numerator  of  the  leg  term  equal  to  zero  produces 
an  expression  for  a  p , ^  ^  of  negative  sign.  When  q«0  this 
is  equivalent  to  p^  of  Eq.  (27).}  Note  also  that  for  q»l, 
as  in  a  case  of  hydrostatic  loading,  all  devlatoric  strains 
are  zero  and  no  p,.  c?n  be  computed.  When  C_  -  0,  p . , 
comes  from  setting  the  numerator  of  the  log  term  equal  to 
zero.  In  that  case  a  valid  p^  1b  obtained  providing 

q  < 

For  the  initial  loading  portion  of  the  devlatoric  curve, 
up  to  s*  (see  Fig.  2),  Eq.  (48)  or  (49)  holds  with  e£°^  ■  0 
and  G  quantities  given  by  Eqs.  (46).  In  the  unloading  stage 

from  s^  to  s^  Eq.  (48)  or  (49)  again  holds  but  e£°^  i« 

*  * 
replaced  by  e^  ,  pq  is  the  pressure  corresponding  to  s^  and 

the  expressions  for  5  are  given  by  Eqs.  (47).  Reloading 

** 

the  stresses  above  s^  may  also  be  computed  from  the 

appropriate  one  of  the  above  two  equations  providing  G 

(o)  ** 

expressions  again  come  from  Eqs.  (46),  and  e^  ■  e^  and 

** 

Po  is  replaced  by  the  pressure  corresponding  to  s^  . 

In  the  special  case  where  pressure  p  exceeds  the 
critical  pressure  of  Eq.  (7)  the  G  quantities  of  Eqs.  (46) 
and  (47)  must  be  redefined  as 
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G  -  G 


G,  - 


i  y_i 

’  *  y  2  +  r 


1-~  /5  Y 

1  +  2q  Y1 


/I 


+  2q 


(qo 


(o) 


„(oK 
°2  > 


g,-o 


for  loading,  while 


-  1  li  4.  ^  ^1U  /  „(<>> 

too  *  GoU  *  4  y2  +  1  +  2q  (q0l 


„(oK 
a2  ) 


L-t  q..  c 

1  +  2q  T1U 


g2  -  0 


for  unloading.  In  either  case,  since  G2  -  0  the  appropriate 
solution  equation  is  Eq.  (49)  for  both  loading  and  unloading 
when  p  >  pc  . 


(51) 


(52) 
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III  MATCHING  LABORATORY  DATA. 

The  true  test  of  any  mathematical  model  is  its  ability 
to  match  experimental  data  for  a  real  material.  The  process 
used  to  choose  the  parameters  such  that  the  model  described 
in  the  previous  section  fits  the  experimental  curves  for  a 
real  soil,  namely  McCormick  Ranch  Sand,  is  discussed  in  the 
present  section.  The  results  for  various  tests  are  compared 
with  typical  experimental  data  in  Section  IV. 

Description  of  the  Available  Experimental  Data. 

The  laboratory  data,  available  at  the  time  the  models 
were  fit,  consisted  of  the  following: 

(a)  A  series  of  static  triaxial  compression  tests  in  which 
stress  difference  was  plotted  against  axial  strain. 

The  tests  were  in  three  groups.  In  the  first  group  the 
load  was  increased  mono tonically  to  failure.  In  the 
second  and  third  groups,  the  specimens  were  loaded  to 
approximately  35Z  and  75Z  of  failure,  respectively,  then 
unloaded  and  reloaded  for  several  cycles  and  finally 
loaded  to  failure. 

(b)  A  composite  static  uniaxial  strain  test:  A  series  of 
both  static  and  dynamic  tests  were  run  in  which  the 
axial  stress  and  axial  strain  were  measured.  The 
radial  stress  was  not  measured.  A  single  uniaxial 
strain  curve  was  estimated  by  W.E.S.  as  the  most 
representative  static  curve  for  the  soil  with  the 


water  content  and  denaity  typical  of  the  triaxial  testa. 

An  Idea  of  the  possible  variation  of  the  uniaxial  stress- 
strain  curve  may  be  seen  in  Fig.  (3)  where  dynamic  curves 
for  three  different  samples,  as  well  as  the  W.E.S.  con¬ 
structed  dynamic  composite  curve,  are  shown. 

(c)  Three  static  hydrostats:  The  actual  measured  values  of 
pressure  and  volumetric  strain  and  the  W.E.S.  constructed 
average  curve  are  shown  in  Fig.  (4).  No  unloading  or 
reloading  was  measured. 

(d)  Three  composite  static  proportional  loading  tests:  The 
tests  were  run  with  the  ratios  of  the  radial  stress  to 
axial  stress  maintained  *t  0.4,  0.6  and  0.8,  respectively. 
The  only  data  available  were  the  three  composite  plots 

of  vertical  stress  versus  vertical  strain.  This  data 
was  used  only  as  a  check. 

(e)  Stress  difference  versus  strain  difference  in  triaxial 
compression.  This  information  was  received  after  the 
original  initial  loading  models  were  already  constructed. 
As  a  result,  for  these  models  the  loading  data  was  only 
used  as  a  check,  however,  the  unloading  data  was  used 

in  constructing  the  unloading  portion  of  the  model. 

The  loading  data  was  used  in  the  improved  models  discussed 
in  Section  IV. 

A  few  comments  are  in  order  to  qualify  the  data.  The 
model  described  in  Section  II  is  assumed  to  apply  to  a  single 
homogeneous  isotropic  (albeit  nonlinear)  material.  Obviously, 
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soil  is  not  homogeneous,  although  in  a  macroscopic  tense  it 

may  behave  as  if  it  were.  However,  different  samples  have 

somewhat  different  properties.  Thus,  at  least  as  tar  as 

initial  loading  is  concerned,  each  test  destroys  the  sample 

so  that  different  samples  must  be  used  in  different  tests. 

It  is  therefore  difficult  in  comparing  results  for  different 

test  configurations  to  decide  whether  the  model  is  in  error, 

or  whether  the  data  for  the  various  tests  simply  apply  to 

different  materials.  An  example  of  this  uncertainty  may  be 

seen  by  comparing  the  initial  bulk  modulus  Kq  =  20  ksi  from 

Fig.  (4)  and  the  initial  constrained  modulus  M  «  K  + 

o  o 

4 

+  —  Go  <  20  ksi  from  Fig.  (3).  The  only  wav  the  two  results 
could  be  compatible  would  be  for  <  0,  an  obvious  absurdity. 

Another  difficulty  was  that  the  data  was  incomplete. 

For  example,  in  the  uniaxial  strain  test  the  lateral  stress 

*) 

Oj  required  to  maintain  zero  lateral  strain  was  not  measured 
If  0^  had  been  measured,  an  independent  pressure-volume 
relation  for  the  sample  in  tht  uniaxial  strain  configuration 
could  have  been  constructed. 

Finally,  it  should  be  noted  that  the  data  came  from  two 
different  laboratories.  The  uniaxial  strain  tests  were  run 
at  W.E.S.  All  other  tests,  the  hydrostats,  the  triaxial  tests 
and  the  proportional  loading  tests,  were  performed  at  Georgia 
Tech  under  W.E.S.  contract.  All  the  Georgia  Tech  samples  were 


At  the  time  the  tests  were  conducted  such  measurements 
could  not  be  made.  Presently,  however,  they  are  made 
as  a  matter  of  course  at  W.E.S. 


prepared  the  same  way  and  tested  on  the  sane  apparatus.  The 
uniaxial  staples  were  prepared  at  W.E.S.  in  a  different 
fashion  and  were  tested  on  a  different  type  of  apparatus. 

A  more  complete  qualification  of  the  data  is  given  by 
J.S.  Zelasko  of  W.E.S.  in  the  preface.  One  more  fact  worth 
mentioning  is  that  the  lateral  deformation  measurements  obtained 
for  specimens  tested  in  the  triaxial  device  were  part  of  an  ex¬ 
perimental  research  and  development  project  Strain  computations 
based  on  these  measurements  were  found  to  be  highly  sensitive  to 
interpretation  of  the  raw  data.  Ref.  [3]. 

Hydrostat-Triaxial  Fit. 

The  problem  of  choosing  material  parameters  to  fit  the 
data  is  greatly  simplified  when  there  are  tests  available  in 
which  as  few  of  tie  independent  variables  as  possible  are 
varied  simultaneously.  The  hydrostat,  in  which  the  pressure 
and  volumetric  strain  are  measured,  is  therefore  a  likely 
candidate  for  computing  the  various  K's  which  appear  in 
Eq.  (15).  A  series  of  pure  shear  tests  in  which  shear  stress 
and  shear  strain  were  measured  and  in  which  the  pressure  was 
kept  constant  during  each  test  would  be  desirable  tests  from 
which  to  determine  the  various  y's  in  the  shear  moduli. 

However,  such  tests  were  not  available. 

Of  the  available  data,  the  stress  difference  at  failure 
in  the  various  triaxial  tests  may  be  used  to  obtain  all  but 
one  constant  in  the  expression  for  GLD  ,  Eqs .  (5)  and  (9). 

At  failure  *  0,  so  that  by  dividing  Eq.  (5)  by  >  0, 
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one  has 


as  a  quadratic'  function  of  pressure  valid  for 


P  <.  P  .  °r 


ao  +  al 


2  Go  .  Y1 


2  2 


p  +  a.p  -  —7-  +  —7-  p  +  —7-  p 


-Yl  -YX 


The  raw  data  of  0,  and  (0,  -  0,)  for  each  test  were 
J  1 _  J  max 


rfLj r=- 


transformed  into  p  an 


polynomial  was  then  fit  through  the  resulting  data. 


decree 


The  results  for  the  various  early  trials  are  shown  in 

Table  I  cases  1  to  9 .  The  second  column  "N"  refers  to  the 

number  of  tests  included  in  the  trial,  while  the  mean 

square  residual  is  a  measure  of  the  width  of  the  scatter 

band.  After  the  three  coefficients  a^  the  ratio  of 

a./a  -  y,/G  and  the  cut-off  pressure  p  =  -a,/2a,  are 
1  o  1  o  c  1  2 


given.  Finally,  the  maximum  va 


rV 

lue  of  V  J2  for 


particular  fit,  the  value  at  p  =  pc  ,  is  given.  Upon 
examination  of  the  output  of  the  first  four  cases  it  was 
evident  that  (a)  systematic  errors  appeared  in  the  fit 
by  trying  to  include  both  high  chamber  pressure  (o^  up  to 
10  ksi)  and  low  chamber  pressure  tests,  and  (b)  there 
appeared  to  be  significant  differences  between  the  noncyclic 


and  the  cyclic  tests, 


It  was  therefore  decided  to  base  the  failure  fit  upon 
the  Low  chamber  pressure  tests  which  were  loaded  straight  to 
failure.  Although  the  mathematical  model,  described  in 
Section  II,  in  which  failure  would  occur  at  the  same  stress 
level  in  both  the  cycled  and  uncycled  triaxial  tests,  ran 


(to  a  certain  extent)  counter  to  the  evidence  in  this  regard, 
no  attempt  was  made  to  alter  the  model  so  that  the  failure 
stress  state  would  depend  upon  the  previous  load  history. 

Sample  16  was  eliminated  from  further  consideration  since 
its  water  content  was  significantly  lower  than  that  of  the 
other  camples  and  its  strength  correspondingly  higher.  Tests 
at  lower  and  lower  chamber  pressures  were  eliminated  (cases 
5-7)  until  all  tests  included  in  the  group  failed  at  pressures 
lower  than  the  computed  p^  for  the  group.  Case  8  which  Included 
thirteen  tests  all  at  chamber  pressures  less  than  or  equal  to 
0.8  ksl  was  this  final  result.  The  pressure  at  maximum  shear 
strength  pc  was  1.166  ksl  and  the  corresponding  value  of 
was  0.2303  ksl.  As  a  citeck,  the  best  constant  through  the 
remaining  uncycled  trlaxlal  data,  case  9,  was  found  to  be 
0.2323  ksl,  which  was  within  lit  of  the  above  value. 

The  trlaxlal  tests  were  also  used  to  obtain  an  estimate 
of  one  of  the  two  "elastic  constants"  Eq  .  This  was  accomplished 
by  extrapolating  backward  to  zero  chamber  pressure  the  best  fit 
of  initial  elope  as  a  function  of  0^  .  The  value  obtained 
using  a  least  squares  procedure  was  Eq  »  14.14  ksl.  The  details 
of  the  procedure  will  be  discussed  later  when  the  second  fit 
based  upon  the  trlaxlal  and  uniaxial  test  results  is  described. 

The  hydrostat,  Fig.  (4),  was  used  in  this  early  fit  to 
obtain  the  values  of  the  various  K's.  The  best  least  squares 
cubic  was  fit  through  the  various  groups  of  data  points.  The 

results  are  shown  in  Table  II.  The  use  of  all  points  slmul- 

2 

taneously,  case  1,  leads  to  nonsense  since  E(y  -  y*)  ,  rather 
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Chan  Che  distance  CO  Che  curve,  is  what  Is  minimized.  Al- 

Chough  there  are  small  differences  in  strain  at  a  given  (high) 

pressure,  Che  large  differences  in  pressure  at  a  given  strain 

drive  the  entire  computation.  Including  the  high  pressures 

even  for  a  single  sample  is  still  inaccurate  for  the  same 

reason.  Terms  higher  than  cubic  should  be  included  If  high 

pressure  data  is  to  be  taken  into  account.  The  last  three 

cases,  each  limited  to  a  single  sample  at  lower  pressures, 

give  meaningful  results.  The  last  two  cases,  in  fact,  are 

very  good  as  shown  by  the  aean  square  residual.  The  quantity 

a  should  be  zero  for  p-0  when  C, .  *  0 .  For  the  last  two 
o  v  kk 

cares  it  is  less  than  2  psi  and  was  neglected.  For  this 
eevly  fit  the  values  of  a^  to  a^  for  case  7  were  used  since 
the  sample  48  data  falls  between  that  of  the  other  two 
samples.  Using  e  -  e^/3  and  the  followin8  values 

were  obtained: 


K  «  a,  »  20.6  ksi 

o  1 

»  6a2  •  -3800  ksi 
K2  -  27a3  -  965000  ksi 

Based  upon  the  value  K  -  20.6  ksi  and  E 

o  o 

obtained  from  the  triaxial  test  initial  slopes 
Poisson's  ratio  and  shear  modulus 


(54) 


14.14  ksi 
the  initial 


v  -0.385 

o 

G  -  5.11  ksi 
o 

are  obtained.  From  G  and  the  coefficients  of  the  triaxial 

o 

failure  fit  given  in  Table  I  the  various  y's  are  computed, 


(55) 


« 
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Yi  -  -Vao  *  '62,9  j 

Y1  "  “Ylal  “  16*05 

y2  “  -Y1a2  "  "  6 . 89/ksi 

The  model  'during  Initial  loading  1b  completely  defined 

by  the  values  given  in  Eqs.  (54)  to  (56).  The  behavior  of 

the  model  in  the  various  test  configurations  was  computed 

*) 

using  specialized  codes  PROP  and  UNAX2  and  compared  with 
the  experimental  data. 


(56) 


There  was  fairly  good  agreement  of  the  computed 
behavior  using  the  present  model  and  the  measured  results 
of  the  triaxial  tests.  This  was  true  of  both  the  stress 
difference  versus  axial  strain  plots  and  the  stress  difference 
varsus  strain  difference  plots,  which  were  received  after  the 
model  (and  the  later  Uniaxial-Trivial  Model  I)  was  already 
constructed.  The  results,  not  shown,  are  generally  com¬ 
parable,  as  far  as  agreement  with  the  triaxial  test  experi¬ 
mental  curves  is  concerned,  with  those  obtained  with  the 
later  Uniaxial-Triaxial  Model  I,  discussed  later.  The  model 
results  are  too  soft,  especially  at  higher  chamber  pressures, 
but  have  the  proper  failure  stress.  Of  course,  the  present 
model,  which  was  based  on  the  hydrostat,  agrees  with  that 
test.  The  results  are  discussed  more  completely  in  Section  IV. 


Unfortunately,  but  not  unexpectedly,  the  present  Hydrostat- 


*) 


Short  descriptions  and  FORTRAN  listings  if  PRO?  and  UNAX2 
are  found  in  the  Appendix. 


Trlaxlal  model  produced  a  stress-strain  curve  in  uniaxial 
strain,  Fig.  (8),  which  was  stiffer  than  the  suggested 
experimental  curve  by  more  than  a  factor  of  two.  This  is 
partially  explained  by  just  the  computed  initial  constrained 
modulus  Mq  -  27.4  ksi  versus  a  value  certainly  less  than 
20  ksi  suggested  by  Fig.  (8). 

Since  the  uniaxial  strain  geometry  is  closer  to  that 
occurring  during  real  events,  and  since  W.E.S.  had  much 
greater  confidence  in  the  uniaxial  data  than  in  the  hydro¬ 
static  data.  It  was  decided  to  forego  the  hydrostat  entirely 
and  base  a  fit  on  the  uniaxial  and  trlaxlal  data  exclusively. 
The  description  of  the  later  model  follows. 

Triaxlal-Uniaxial  Fit. 

The  present  approach  utilized  the  failure  envelope  from 

the  trlaxlal  tests  to  obtain  the  ratio  of  the  y's,  the 

initial  slopes  of  the  trlaxlal  tests  to  estimate  E  and 

o 

the  initial  slope  of  the  uniaxial  test  M  as  the  other 

o 

"elastic"  constant.  Finally,  and  Kj  were  found  by  a 
trial  and  error  routine  to  obtain  the  properly  shaped  uni¬ 
axial  strain  curve.  This  later  approach  was  done  In  a  more 
systematic  fashion  since  experience  had  already  been  gained 
during  the  earlier  fit. 

Moreover,  the  data  was  somewhat  refined  from  what  was 
used  previously.  For  example,  whereas  before  one  had  to 
select  which  trlaxlal  failure  data  to  use,  presently  all  the 
test  samples  which  had  water  contents  and  chamber  pressures 


In  the  Appropriate  range  were  used.  In  addition,  the 
numerical  values  were  slightly  sore  accurate.  The  data 
used  to  compute  the  new  failure  envelope  is  shown  in 


Table  III.  The  raw  data.  0,  and  (0,  -  0,)  ,  which  came 

’3  13  max  ’ 


The  uppermost 


from  W.E.S.  were  used  to  compute  p  and 
section  of  the  table  are  the  uncycled  tests:  the  middle  and 
bottom  sections  are  the  tests  in  which  unloading  took  place 
at  35%  and  75%  of  the  failure  stress,  respectively. 


The  results  of  the  various  least  squares  fits  are  given 
in  Table  I,  cases  10  to  13.  The  data  from  the  uncycled  tests, 
case  10,  contain  much  less  scatter  than  the  tests  with  un¬ 
loading  as  can  be  seen  from  the  mean  square  residual.  This 
is  also  apparent  from  Figs.  (5)  and  (6).  In  Fig.  (5)  the 


failure  fit  {  R  versus  p)  is  plotted  for  the  uncycled 
tests  together  with  the  data  points.  The  scatter  band  is 
quite  small)  Fig.  (6),  on  the  other  hand,  which  is  the  same 
plot  for  each  of  the  cyclic  loading  sets,  shows  a  large 
scatter  band,  particularly  at  higher  pressure.  It  is  noted, 
however,  that  the  maximum  value  of  r  is  approximately  the 
same  for  all  the  groups  of  tests.  The  model  was  based  on 
case  10,  the  uncycled  triaxial  tests  because  of  the  small 
scatter  band.  It  is  seen  that  the  coefficients  a^  did  not 
change  drastically  from  the  previous  model,  case  8. 


As  in  the  earlier  Hydros tat-Triaxial  fit,  one  of  the 
"elastic  constants"  was  obtained  from  the  initial  slopes  of 
the  triaxial  tests.  The  results  are  summarized  in  Table  IV. 
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Various  degree  polynomial  f'its  were  tried  on  various  segments 

of  the  data.  The  emphasis  was  placed  on  the  lower  pressure 

(o j  <_  400  psi)  tests  since  the  primary  interest  was  an 

accurate  estimate  of  E  *  a  .  There  was  less  scatter  in 

o  o 

the  cyclic  test  data  than  in  the  noncyclic  tests,  as 

witnessed  by  comparing  the  mean  square  residual  of  cases  5 

and  6  with  case  4.  However,  since  there  was  no  reason  why 

the  Initial  slope,  prior  to  any  unloading,  should  depend  on 

the  subsequent  loading  history,  all  tests  were  included. 

The  result  used,  case  10,  a  linear  fit  through  all  the  tests 

with  O,  <  400  psi,  was  E  -  12.2  ksi.  The  last  column  in 
J  —  o 

Table  TV  shows  the  computed  value  of  E  at  “  100  psi.  The 
computed  value  for  case  10,  30  ksi,  falls  within  the  raw  data. 

The  measured  values  of  initial  slope  together  with  the  fit 
used  are  plotted  against  chamber  pressure  in  Fig.  (7).  The 
computed  curve  is  dashed  above  0.4  ksi  since  it  is  not 
applicable.  The  wide  scatter  in  the  initial  slope  data  is 
immediately  apparent  from  Fig.  (7). 

The  remaining  "elastic"  constant  was  the  initial  slope 

of  the  uniaxial  strain  test  M  .  This  was  taken  as  16.5  ksi 

o 

based  on  the  experimental  curve.  Fig.  (8),  suggested  by  W.E.S. 
Using  the  relations  between  the  elastic  constants 

M  -  K  +  G  (57) 

and 

9KG 


E 


3K  +  G 


(58) 


one  may  solve  for  G  and  K  In  terns  of  E  and  M  as 


G 


(3M  +  E)  -  V  (9M  -  E) (M  -  E) 
8 


(59) 


and 


K  -  M  -  G  (60) 

With  E  «  12.2  ksi  and  M  «  16.5  ksit  Eqs.  ( )  and  (60) 
o  o 

give 


G  -  4.69  ksi 

o 

K  -  10.24  ksi 
o 


(61) 


The  Initial  value  of  Poisson's  ratio  V  turned  out  to  be  0.30. 

o 

From  the  values  of  a^  of  the  failure  envelope  (case  10 
In  Table  I)  and  Gq  the  various  y's  were  computed. 


-G  /a  « 

-64.2 

0  0 

-Vi  - 

18.9 

-*la2  “ 

-8. 76/ksi 

(62) 


The  remaining  two  constants  required  to  complete  the 
material  description  for  Initial  loading,  and  Kj  ,  were 
obtained  by  a  trial  and  error  scheme  which  was  continued 
until  a  suitable  fit  to  the  experimental  uniaxial  test  was 
obtained.  First,  with  some  trial  value  of  held  constant, 
was  varied  so  that  the  Influence  of  a  change  in  on  the 
stress-strain  curve  could  be  ascertained.  Then  was 
adjusted  to  make  the  computed  curve  agree  with  the  experi¬ 
mental  one  in  the  middle  stress  range  (o^  between  300  and 
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500  psi).  Next,  with  held  constant,  K2  was  varied  and 
adjusted  to  make  the  curves  coincide  in  the  high  stress 
range.  With  this  value  of  ,  the  cycle  was  repeated  to 
"fine  tune"  .  The  final  values 


K:  =  -1,250  ksi 
K2  =  97,000  ksi 


(63) 


were  obtained  in  a  few  Iterations.  The  computed  curve  is 
superimposed  on  the  experimental  curve  in  Fig.  (8). 


After  all  the  material  parameters  for  initial  loading 
were  determined,  plots  of  stress  difference  versus  strain 
difference  for  some  of  the  triaxial  tests  ware  received. 

These  plots  would  have  enabled  one  to  obtain  G  ,  Y,  and  y, 
directly  from  the  initial  slope  as  a  function  of  .  The 
resulting  shear  modulus  would  have  been  significantly  stiffer 
than  the  one  obtained  from  the  failure  envelope  and  E 

o 

This  is  shown  clearly  in  Table  V  where  the  computed  values 

labeled  Uniax-Triax  I  are  all  smaller  than  the  lower  bounds 

of  the  measured  values.  (The  computed  values  of  2G.  , 

r  initial 

for  the  improved  models  discussed  in  Section  IV,  also  given 
in  Table  V,  agree  much  better  with  the  measured  values.) 


It  is  not  clear  whether  this  discrepancy  is  in  the 
model  formulation  or  in  the  measurement  of  lateral  strain. 
In  any  v-.-.se,  this  discrepancy  shows  up  in  the  results  which 
will  be  discussed  later. 


Although  the  plots  of  stress  difference  versus  strain 
difference  were  not  used  at  all  in  determining  the  present 


29 


* 


*) 

loading  model  ,  they  were  the  basis  upon  which  the  unloading 
shear  modulus  parameters  were  determined. 


Unloading. 

fhe  values  of  the  various  material  parameters  which 
describe  the  model  in  unloading  were  determined  completely 
from  the  cycled  triaxial  tests.  Where  there  were  several 
unload-reload  cycles,  the  initial  unloading  segment  was  used. 
Since  values  of  the  shear  modulus  were  directly  obtainable 
from  the  slopes  of  the  versus  curves,  which 

had  just  become  available,  these  slopes  were  the  basis  of  the 
parameters  in  G^  .  After  was  determined,  the  measured 

slopes  of  the  conventional  triaxial  plots  E  were  used  to 
calculate  .  The  unloading  portion  of  the  uniaxial 

strain  test  was  used  only  as  a  check. 


The  slope  of  the  stress  difference  versus  strain 

difference  plots  were  measured  at  the  top  and  bottom  of  the 

initial  unloading  portion  of  each  of  the  cyclic  triaxial  tests 

*  *  **  ** 

available,  i.e.,  at  the  points  (s^  ,  e^)  and  (s ^  ,  e^  ), 

respectively,  in  Fig.  (2).  The  stress  difference  (cr^  -  0^) 

at  each  end  and  the  chamber  pressure  0^  were  noted.  From 

(o^  -  <Jj)  and  0 ^  values  of  p  and  ~\J J 2  were  computed.  The 

results  are  summarized  for  the  tests  unloaded  from  75 7.  and 

35%  of  failure  in  Tables  VI  and  VII,  respectively.  At  the 

bottom  the  dominant  effect  is  due  to  p  since  yr  is  small. 

The  most  self-consistent  values  of  2G„  occur  in  the 

Bottom 


*) 


They  were  uted  in  formulating  the  improved  models 
Uniax-Triax  II  and  Hydro-Triax  II  discussed  in  Section  IV. 
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samples  unloaded  from  75%  of  failure  with  0^  of  200  and 
400  psi.  These  values  are  boxed-in  in  Table  VI.  By 
assuming  V J ?  is  negligible  and  using  Eq .  (8)  to  relate 
Y1r  and  Y?„  ,  the  expression  for  GU(J  ,  Eq .  (6),  is  reduced 


to  one  involving  two  independent  parameters  G  and  Y...  . 

r  oU  ' 1U 

Using  the  average  value  of  p  and  the  average  value  of 
^GUX  ^or  tl'e  two  k°xe£*  cases,  one  has  two  simultaneous 


equations  for  G  and  Y^  • 


At  the  top  of  the  unloading  portion 


f. r 


is  clearly 


important.  The  most  self-consistent  values  of  2G„  occur 
K  Top 

in  the  352  group  with  “  100  psi.  Only  samples  112  and 

115  are  boxed— in  in  Table  VII,  since  sample  108  obviously 

belongs  in  the  752  cyclic  group  even  though  it  is  labeled 

352.  Using  the  average  2GTop  of  samples  112  and  115,  and 

the  known  values  of  p  ana  together  with  the  previously 

determined  values  GqU  and  Y ^ u  (and  ^2^)  leads  to  the 

determination  of  y,,,  • 

1  u 


The  tentative  values  of  the  unloading  shear  parameters 
so  determined  are 


G 

Y 

Y 

Y 


oU 

1U 

2U 
i  U 


7.0  ks  i 
51.1 

-23. 7/ksi 
172 


(64) 


The  values  of  the  unloading  shear  modulus  computed  with  the 
constants,  Eqs .  (64),  (slightly  rounded  off)  are  compared 
with  the  range  and  average  measured  values  for  each  case 


in  Tables  VIII  and  IX.  Fairly  good  agreement  is  seen  for  2Gtlu 
at  the  bottom  of  the  unloading  segment,  Table  VIII,  with  only 
three  cases  falling  slightly  outside  the  scatterband.  At  the 
top  of  the  unloading  segment.  Table  IX,  however,  agreement  is 
found  only  for  a  single  case,  namely,  the  35X  -  100  psi  group 
used  to  compute  y^  in  the  first  place.  All  other  computed 
values  of  2GyN  fall  significantly  below  the  bottom  of  the 
scatterband,  indicating  a  choice  of  y^y  which  was  far  too  low. 


New  higher  values  of  y^y  were  used  to  compute  new 
calculated  values  in  Tables  VIII  and  IX.  Better  agreement 
was  obtained  at  the  top,  but  the  higher  values  of  y  also 
caused  higher  calculated  values  of  2Gyjj  at  the  bottom.  As  a 
result,  slight  modifications  were  made  in  Goy  and  y  (and  y^y 
to  improve  the  agreement  at  the  bottom.  The  final  values  of 
the  unloading  shear  parameters  used 


G  „  »  6.0  ksi 

oU 

TlU  "  40’° 

Y2u  “  -18.5/ksi 

*iu  '  500 


\ 


J 


represent  a  reasonable  attempt  to  stay  within  the  scatter. 
The  choice  is  not  claimed  to  be  unique.  The  values  of  2GyN 
computed  with  the  parameters,  Eq3.  (65),  are  compared  with 
the  measured  values  in  Tables  VIII  and  IX.  At  the  bottom. 
Table  VIII,  three  cases  fall  outside  the  range  of  measured 
values ,  the  same  as  for  the  tentative  values,  Eqs.  (64). 
However,  at  the  top  of  the  unloading  segment,  Table  IX,  the 
calculated  values  fall  outside  the  range  of  measured  values 
for  only  two  cases  (including  the  35%  -  100  psi  case  which 
agreed  previously!). 


(65) 


In  order  to  compute  the  parameters  KqU  and  In  the 

unloading  bulk  modulus  the  slopes  of  the  conventional 
triaxial  plots  E  were  used  in  conjunction  with  the  valuer 
Eqs.  165).  The  slopes  d(a^  -  O^/de^  »  E  were  measured  near 
the  bottom  of  the  initial  unloading  segment  of  representative 
samples.  The  results  are  summarized  in  Table  X.  From  the 
chamber  pressure  O and  the  value  of  (0^  -  a^)  at  which  E^N 
is  measured,  p  and  fZ  may  be  computed,  and  with  the  help 
of  values,  Eqs.  (65),  G^,  as  well,  Firally,  from  Eq .  (58) 


EC 


9G  -  3E 


(66) 


so  that  may  be  computed  for  each  sample  listed  in 

Table  X.  It  is  apparent  that  the  two  values  of  at 

p  «  0.112  ksi  are  in.. -nsistent ,  whereas  the  pair  of  values 
of  at  p  -  0.22  ksi  and  at  p  =  0.42  ksi  agree  quite  well. 

Taking  the  average  pressure  and  corresponding  values  of  K 
for  those  two  cases,  the  two  constants  appearing  in  Eq.  (4), 
and  ,  are  found  by  solving  two  simultaneous 

equations  to  be 


K 


oU 


K 


1U 


32.0  ksi 
143 


(67) 


The  values  Eqs.  (65/  and  (67)  completely  define  the 
model  in  unloading.  It  should  be  noted  .'.hat  the  in¬ 
equalities  Eqs.  (18),  (19),  (21),  (22)  and  (24)  are  satisfied 
by  the  unloading  values  Eqs.  (65)  and  (67)  and  either  the 
Hydros tat-Triaxial  fit  [Eqs.  (54),  (55)  and  (56)]  or  the 
I'ni  a  X  ia  1-Tr  i  axi  a  1  fit  [Eqs.  (61),  (62)  and  (63)]. 
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IV  DISCUSSION  OF  RESULTS. 

Initial  Loading. 

The  behavior  of  the  model,  in  particular  the  Uniaxial- 
Trlaxial  fit,  in  initial  loading  is  illustrated  in  Figs. 

(8) -(20).  The  stress-strain  relation  in  uniaxial  strain, 

Fig.  (8),  closely  follows  the  typical  experimental  curve 
suggested  by  U.E.S.  The  maximum  deviation  is  less  than 
15  psi  for  stresses  less  than  500  psi  and  about  25  psi  at 

0^  -  700  psi.  This  is  not  surprising  since  the  experimental 
curve  was  not  used  to  construct  the  model.  Also  shown  in 
Fig.  (8)  is  the  stress-strain  curve  predicted  by  the 
earlier  Hydrostat-Triaxial  fit.  As  mentioned  previously, 
the  Hydrostat-Triaxial  fit  results  in  a  uniaxial  strain 
stress-strain  relation  which  is  far  stiffer  than  any  observed. 

The  stress  difference  versus  axial  strain  results  in 
triaxial  compression  tests,  loaded  to  failure,  are  given 
in  Figs.  (9)-(12)  for  chamber  pressures  of  0^  -  0.1,  0.2, 

0.4  and  0.8  ksi.  The  plots  of  stress  difference  versus 
strain  difference  for  the  same  configuration  and  the  same 
chamber  pressures  are  shown  in  Figs.  (13)-(16).  The  various 
experimental  curves,  several  for  each  figure,  are  the  light 
8oxji.d  lines.  All  the  available  experimental  curves  are 
included,  the  cycled  tests  as  well  as  the  uncycled  ones. 

The  unloading-reloading  cycles  have  been  left  out  in  Figs. 

(9) -(12)  for  clarity.  The  single  cycle  of  unloading-reloading 
vas  included  in  Figs.  (13)-(15)  to  give  some  idea  of  the 
spread  in  cyclic  data.  No  such  data  is  shown  in  Fig.  (16), 
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again  for  clarity.  It  should  be  noted  that  the  exact  values 

of  initial  slopes  and  (to  a  much  lesser  extent)  stress 

differences  at  failure  differ  from  those  given  in  Fig.  (7) 

and  Table  III.  The  actual  fit  was  constructed  from  uncorrected 

raw  test  data,  while  the  experimental  curves  shown  represent 

processed  data.  The  scatter  in  the  initial  modulus  E.  .  .  , 

initial 

is  reduced  significantly  in  Figs.  (9)-(12)  as  compared  with 
Fig.  (7). 

The  heavy  solid  curve,  labeled  "Uniax-Triax  1",  is  the 
curve  computed  using  the  Uniaxial-Triaxial  fit  discussed  in 
Section  III.  The  earlier  Hydros tat-Triaxial  fit,  not  shown, 
gave  results  which  were  very  similar.  The  computed  stress 
difference  at  failure  is  shown  as  a  horizontal  asymptote. 

The  other  heavy  curves  labeled  "Uniax-Triax  2"  and  "Hydro- 
Triax  2"  are  Improved  fits  which  are  discussed  below. 

With  the  exception  of  Figs.  (9)  and  (13)  which  correspond 

to  the  tests  run  at  0^  =  100  psi,  all  computed  curves  approach 

the  correct  stress  at  failure.  The  reason  for  the  computed 

(0,  -  0_)  being  high  at  0,  ■  100  psi  is  evident  from 
1  J  max  j 

Figs.  (5)  and  (6)  where  all  the  data  points  at  100  psi  fall 
below  the  failure  fit.  Although,  in  general,  the  computed 
curves  on  the  Uniax-Triax  1  fit  approach  the  proper  value  at 
failure,  they  are  too  soft  at  stresses  significantly  below 
failure.  The  discrepancy  appears  to  increase  with  chamber 
pressure  in  both  the  conventional  triaxial  and  stress 
dif ference  -  s train  difference  plots.  The  similar  trend  in 
both  sets  of  plots  suggests  that  the  shear  modulus  GLI)(p, 
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is  at  fault  and  that  altering  the  values  of  the  parameters 
Gq  ,  Yi  .  Y2  Yj  may  lead  to  significant  improvements. 

The  proportional  loading  teats  are  shown  in  Figs.  (17), 

(18)  and  (19)  for  stress  ratios  q  of  0.4,  0.6  and  0.8, 

respectively.  Since  these  tests  were  in  no  way  used  to 

construct  the  model,  they  serve  as  a  check  on  its  validity 

in  a  different  loading  geometry.  For  q  -  0.4,  Fig.  (17), 

the  computed  axial  stress  at  failure  agrees  well  with  the 

peak  stress  in  the  experimental  curves.  As  q  Increases, 

the  computed  peak  axial  stress  becomes  larger  than  that 

which  is  measured.  At  a  stress  ratio  of  q  *  0.8,  Fig.  (19), 

the  discrepancy  is  quite  pronounced,  the  computed  value  of 

0^  being  2.01  ksl  while  the  measured  peak  stress  ranges 
max 

from  1.22  to  1.62  ksl. 


It  turns  out  that  the  axial  stress  at  failure  is  quite 
sensitive  to  uncertainties  in  the  stress  ratio  q.  Since 
these  stress  levels  correspond  to  the  horizontal  (or  nearly 
so)  portion  of  the  failure  surface.  Fig.  (5), 


_  (CT1  ~  q3)max  _  0,403 

1  1  -  q  1  -  q 

max 


ks  1 


If  when  the  test  was  run  q  had  been  0.75  rather  than  0.80, 

(a  not  unreasonable  error  of  6%),  according  to  Eq.  (68) 

0^  would  have  been  1.61  ksl,  a  value  very  close  to  one 
max 

actually  measured.  Even  the  lowest  value  0^  «  1.22  ksi. 


max 

corresponds  to  a  value  q  0.67,  only  16Z  below  the  nominal 

value.  Of  course,  the  present  discussion  applies  to  a  fixed 

value  (0,  -  0.)  .  As  an  example  of  the  possible  scatter 

1  3  max  r  r 


(68) 
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from  sample  to  sample  in  the  stress  difference  at  failure 
at  stress  levels  of  =  1.2  ksi  refer  to  Fig.  (12).  The 
shape  of  the  curves,  labeled  "Uniax-Triax  1"  in  Figs. 

(1  7)— (19)  agrees  reasonably  well  with  those  measured. 

Finally,  the  hydrostatic  test  results  are  shown  in 
Fig.  (20).  The  curve  marked  "Hydro-Tr iax  Fit"  agrees  well 
with  the  data,  at  least  up  to  800  psi.  This  is  no  surprise 
since  it  was  constructed  to  do  so,  see  Section  III.  The 
curves  based  upon  the  uniaxial  test,  on  the  other  hand, 
both  the  present  Ur.iax-Triax  1  fit  and  the  improved  Uniax- 
Triax  2  discussed  below,  are  considerably  softer  than  the 
measured  results. 

Un  loading- Re loading. 

The  behavior  of  the  model  in  unloading  and  reloading 
is  illustrated  in  Fig.  (8)  in  the  uniaxial  strain  test  and 
in  Figs.  (21)-(28)  in  the  triaxial  configuration.  In  the 
uniaxial  test  only  unloading  was  measured,  while  in  the 
triaxial  case  complete  cyclic  data  was  available. 

The  computed  unloading  curve  in  uniaxial  strain,  Fig. 
(8),  agrees  reasonably  well  with  the  measured  curve,  except 
for  the  sharp  tail  at  the  bottom.  Of  course,  the  scale, 
chosen  to  portray  the  loading  portion  well,  more  or  less 
obscures  small  differences  in  the  nearly  vertical  unloading 
portion.  The  tail  at  the  bottom  of  the  unloading  stroke  has 
been  impossible  to  reproduce  exactly  with  any  of  a  variety 
of  mathematical  models  with  which  the  author  has  had 
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experience.  Upon  reloading,  a  hysteresis  loop,  albeit  a 
small  one  is  formed  and  the  reloading  portion  then  crosses 
the  unloading  portion  and  approaches  the  extended  initial 
loading  curve. 

The  results  of  the  cyclic  triaxial  tests,  unloaded  from 
75%  of  the  failure  load  are  shown  in  Figs.  (21)-(24);  the 
35%  cyclic  tests  are  shown  in  Figs.  (25)-(28).  The  con¬ 
ventional  plots,  stress  difference  versus  axial  strain  are 
in  part  (a),  while  the  stress  difference  versus  strain 
difference  plots  are  in  part  (b) .  To  avoid  confusion,  only 
a  single  experimental  curve  is  shown  in  each  figure. 

The  curve  chosen,  in  each  case,  is  typical  of  the 
better  data.  It  should  be  noted  that  the  starting  point 
for  unloading  is  determined  by  the  (initial)  loading  model 
and  not  the  unloading  fit. 

The  unloading  behavior  in  the  triaxial  te:t  agrees 
reasonably  well  with  the  experimental  results,  especially 
when  the  wide  scatter,  c.f.  Tables  VI  and  VII,  in  the  latter 
is  considered.  On  reloading,  however,  the  model  is  clearly 
too  soft.  Upon  multiple  cycles  of  J oadlng-unloading,  the 
stress  difference  versus  strain  difference  plots  predicted 
by  Eqs .  (5)  and  (6)  [or  Eqs.  (9)  and  (10)  where  applicable] 
would  consist  of  multiple  reproductions  of  the  first  cycle, 
each  displaced  in  the  strain  direction.  In  the  corresponding 
stress  difference  versus  axial  strain  plots,  some  difference 
would  exist  betvaen  the  first  (IniMal)  loading  and  subsequent 
reloading  since,  while  G  «  G  ,  K  «  K  during  reloading. 

Li  V  UN 
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Nevertheless,  the  model  predicts  much  more  "walking”,  i.e., 
excessive  strain,  than  is  observed. 

Suggested  Modifications. 

The  I’niax-T.'iax  Fit  1  model  described  above  adequately 
predicts  the  behavior  of  McCormick  Ranch  Sand  in  the  various 
available  laboratory  tests  with  three  exceptions.  (1)  The 
hydrostat  is  too  soft;  (2)  the  initial  loading  behavior  in 
triaxial  tests  is  somewhat  too  soft;  and  (3)  on  reloading 
the  model  predicts  excessive  "walking".  As  far  as  (1)  is 
concerned,  the  experimental  data  is  simply  inconsistent  and 
no  modification  of  the  existing  theory  can  eliminate  the 
discrepancy  between  the  measured  uniaxial  strain  and  hydro¬ 
static  behavior.  The  remaining  discrepancies  can  be  largely 
eliminated  by  the  modification  suggested  below. 

The  general  trend  is  for  the  bulk  modulus  to  Increase 
with  increasing  stress  level  (beyond  some  low  stress)  and 
for  the  shear  modulus  to  increase  with  pressure  and  decrease 
with  increasing  stress  deviators.  In  a  strain  controlled 
test,  such  as  the  uniaxial  strain  test,  the  pressure  volume 
relation  dominates,  since  the  devlatoric  strains  are  con¬ 
strained  by  geometry.  The  bulk  modulus  increases  much 
faster  than  the  shear  modulus  so  that  K/G  becomes  much 
larger  than  unity.  The  slope,  or  constrained  modulus 

M  »  K(1  +  3  |)  =  K 

so  that  at  moderately  high  stresses  the  stress-strain 


(69) 


relation 


is 


insensitive  to  G. 


Ia  a  stresa  controlled  teat,  such  as  the  triaxial 
compreaalon  test,  where  the  deviatorlc  stress  increases 
faster  than  the  pressure,  the  ratio  G/K  becomes  small  since 
G  approaches  zero.  The  slope  at  higher  stress  differences 


E 


3C 


is  thus  insensitive  to  K.  Of  course,  the  slope  of  the 
stress  difference  versus  strain  difference  curves  are 


exactly  2G  and  are  independent  of  K.  The  slope  do^/ds^ 
in  a  proportional  loading  test  is  intermediate  between 
Eqs .  (69)  and  (70) 


(70) 


d£l  3( 1-q)  +  |  (1  +  2q ) 


and  depends  on  the  stress  ratio  q.  The  bulk  modulus 
becomes  more  Important  as  q  approaches  one. 


(71) 


In  view  of  the  discussion  above,  the  fact  that  the 
(Of  -  o3>  versus  both  and  (e^  -  triaxial  test  plots 

are  too  soft  clearly  suggests  that  the  values  of  Gq  ,  y^  , 

Y^  and  y ^  given  by  Eqs.  (61)  and  (62)  are  in  error. 

In  Fig.  (29),  the  measured  values  of  initial  slopes 
2GLD  are  plotted  against  chamber  pressure.  A  fit  through 
the  data,  also  shown,  intersects  the  2G  axis  at  approximately 
lb  ksi.  Using  Gq  *=  8.0  ksi  and  the  original  failure  fit 
(Case  10,  Table  I),  the  modified  shear  modulus  parameters 


become 
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G  -  8.0  ksl 

o 

Yl  -  -110.0 
y1  -  32.4 

f 2  m  ”15.0 /ksl 


(72) 


The  function  2G^p(p,C)  based  on  Eqs.  (72)  Is  shown  In 
Fig.  (29),  as  is  that  based  on  Eqs.  (61)  and  (62).  The 
stlffer  shear  modulus  cones  much  closer  to  the  measured 
values . 


Two  modified  complete  initial  loading  models  were  used 
to  recompute  the  various  tests  for  which  measurements  were 
available,  Figs  (9)  to  (20)j  both  used  Eqs.  (72)  to  define 
the  shear  modulus.  The  "Hydro-Triax  Fit  2"  used  the  bulk 
constants,  Eqs.  (54),  found  previously  based  on  the  hydro¬ 
static  tost.  The  second  model,  the  "Unlax-Triax  Fit  2", 

used  the  initial  constrained  modulus  M  *  16.5  ksl  and 

o 

repeated  the  trial  and  error  procedure  used  previously  to 
fir„  and  so  that  the  uniaxial  strain  test  was 
adequately  reproduced.  The  final  values  of  the  bulk 
constants  were 


K  »  5 . 8  3  ks  i 
o 

K  -  80  ksi 
K2  -  30,000  ksi 

The  computed  uniaxial  strain  curve  for  the  Uniax-Triax  Fit  2 
is  compared  with  the  experimental  curve  in  Fig.  (30).  The 
agreement  is  at  least  as  good  as  that  in  Fig.  (8),  the 
original  uniaxial-triaxial  fit.  Also  shown  are  the  uniaxial 


(73) 
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curves  produced  by  both  Hydros tat-Triaxial  fits.  The 
closeness  of  the  two  clearly  validates  the  contention, 

Eq.  (69),  that  the  uniaxial  strain  curve  Is  Insensitive 
to  the  shear  modulus.  The  unloading  model  used  to  compute 
the  unloading-reloading  portion  of  Fig.  (30)  was  that  given 
earlier  by  Eqs.  (65)  and  (67)  with  the  exception  that  GqU 
was  increased 

G  -  8.0  ksi 
o  u 

so  that  inequality,  Eq.  (21),  lj  not  violated. 

It  is  worth  noting  that,  although  the  two  Uniaxial- 
Trlaxial  fits  give  substantially  the  same  stress-strain 
curve  in  uniaxial  strain,  the  corresponding  stress  paths 
differ  appreciably,  see  Fig.  (31).  In  fact,  the  initial 
Poisson's  ratio  in  the  Uniix-Triax  Fit  2  is  only  0.029 
(perhaps  unrealistically  low)  although  by  an  axial  stress 
of  100  psi  it  is  already  0.23.  The  radial  stress  was  not 
measured  so  /:hat  the  actual  stress  path  remains  unknown. 

The  stress  difference  versus  axial  strain  curves  in 
tiiaxial  compression  produced  by  both  modified  fits  are 
shown  in  Figs.  (9)-(12).  As  suggested  by  Eq.  (70),  the 
curves  are  close  together  and  thus  insensitive  to  the  bulk 
relation.  Of  course,  the  curves  in  Figs.  (13)-(16)  depend 
only  on  G  and  the  two  modified  fits  are  identical.  It  is 
seen  from  Figs.  (9)-(16)  that  the  shear  modulus  parameters, 
Eqs.  (72)  produce  results  which  are  in  excellent  agreement 
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with  the  data. 


In  the  proportional  loading  tests.  Figs.  (17)— (19) , 
the  two  modified  fits  grow  further  apart  as  q  increases,  as 
suggested  by  Eq .  (71).  The  Hydro-Triax  Fit  2  becomes  much 
too  stiff  when  q  «  0.8.  The  Untax-Triax  Fit  2  agrees  slightly 
less  well  with  the  data  than  the  softer  Unlax-Triax  Fit  1; 
nevertheless,  the  agreement  with  the  data  is  not  unreasonable. 
Finally,  the  hydrostats  based  on  all  fits  are  shown  in 
Fig.  (20).  Both  uniaxial  fits  are  much  softer  than  the 
measured  results,  a  discrepancy  discussed  earlier. 


The  present  model,  which  makes  no  distinction  between 
initial  loading  in  shear  and  subsequent  reloading  in  shear, 
results  in  excessive  strains  di’eing  cyclic  loading.  One 
possibility,  analogous  to  the  bulk  modulus  treatment,  woulu 
be  to  use 


.  * 

when  J0  >  0 


■ 

and 


J 


2 

max 


P  ’ 


.i  it 

when  <  0  or  <  Jj 

max 


t  i 

where  is  the  maximum  previous  value  of  •  Use  of 

max 

Eqs.  (75)  would  eliminate  all  additional  hysteresis  and 
additional  permanent  strain  after  the  first  cycle  of  cyclic 
loading.  Another  possibility  is  the  introduction  of  a  new 
function  G__  which  would  be  a  linear  combination  of  the  two 

K  Ei 

independent  shear  moduli  G^n  and  G^  ,  namely 


(75) 


where  a,  8  and  n  are  constants,  and  where  F/F  is  a  measure 

111 

of  the  present  state  of  stress  relative  to  the  maximum  (in 

some  sense)  previous  state  of  stress.  When  F/F  -  0,  for  the 

in 

i 

behavior  to  be  continuous  when  J,  (or  F)  »  0,  G__  ■  G„u  , 

4  RL  UN 

while  when  F/F  »  1  one  sets 
m 

°RE  *  “  0l»  +  (1-6>  °U»  *  °LD  <77) 

This  produces  the  relation  between  a  and  8 

(1-“)  WV  -  (1*B)  CU*(V  <  8> 

which  varies  with  the  loading  history  of  each  point.  How¬ 
ever,  Eq.  (78)  is  satisfied  identically  for  all  loading 
histories  if  a  *  8  -  1. 


The  measure  of  the  state  of  stress  will  be  related  to 


GLD<P’J2>  >  ° 

CLD(p,J2) 

gld(p.0) 

I  I 

When  J2  ■  0,  F  =  0,  while  if  J2  is  such  that  G^p  ■  0,  F  ■  1, 
The  function  F  is  thus  restricted  to  the  range 


(79) 


0  <  F  <  1  (80) 

The  final  form  of  GOI.  is 

CRE  '  <f'“  °LD<’-V  *  11  -  ‘f7”1  W-V  (81> 

m  m 

with  F  given  by  Eq.  (79)  and  F^  its  maximum  previous  value 
at  the  point. 
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The  reloading  curve  will  atart  out  with  the  unloading 
slope.  For  small  (relative  to  the  previous  maximum)  cyclic 
loading  the  hysteresis  and  "walking"  will  be  small.  Since 
Gl'X  -  GLD  •  CRE  i*  restricted  to  the  range 


CLD  ±  C*E  i  CUN 


(82) 


As  the  load  is  increased  towards  its  maximum  previous  value, 
there  is  a  smooth  transition  to  the  virgin  loading  curve. 

The  suddenness  of  the  transition  is  controlled  by  the 
parameter  n;  for  F/F  <  1  as  n  ■»  •  ,  G_„  •+■  G„„  ,  while  as 
n  0,  G  -*■  G.  n  for  F  >  0.  The  value  of  n  could  be  chosen 
so  that  the  reloading  shear  stress  -  snear  strain  curve  is 
always  concave  downward,  thereby  eliminating  the  possibility 
of  any  shear  shocks.  Work  on  reloading  is  continuing. 
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v  CONCLUSIONS  AND  RECOMMENDATIONS. 

A  model  has  been  developed,  naaely  the  Unlax-Trlax  Fit  2, 
which  agrees  excellently  with  the  available  Initial  loading 
uniaxial  strain  and  trlaxlal  compression  data.  The  agreement 
with  the  proportional  loading  data  and  the  unloading  trlaxlal 
data  Is  adequate.  The  questions  of  the  hydrostat  and  reloading 
In  shear  are  discussed  in  Section  IV. 

On  the  basis  of  the  present  study,  the  essential  In¬ 
gredients  for  fitting  a  variable  aodulus  nodel  of  the  type 
described  In  Section  II  to  a  given  noil  are  (a)  the  stress 
strain  relation  In  uniaxial  strain,  (b)  the  failure  envelope 
and  (c)  the  stress  difference  -  strain  difference  relation 
in  a  teat  where  shear  effects  predominate,  such  as  a  trlaxlal 
compression  test.  The  measurement  of  the  lateral  stress  in 
a  uniaxial  strain  test  Is  highly  desirable.  Also  desirable 
Is  a  hydrostatic  ter t  which  contains  unloading-reloading  and 
which  Is  compatible  with  the  uniaxial  results. 

The  question  of  reloading  in  shear,  e.g.,  finding  a 
reasonable  value  for  the  parameter  n  in  Eq.  (81),  is  presently 
left  unanswered. 
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TABLE  III 


McCORMICK.  RANCH  SAND  -  TRIAXIAL  TEST  FAILURE  DATA 


— 

SAMPLE 

°3 

°1  ’  °3 

P 

t 

J2 

88 

0*100 

0.202 

0.167 

0.116 

76 

0.100 

0.191 

0.163 

0.110 

91 

0.100 

0.198 

0.166 

0.114 

97 

0.100 

0.197 

0.165 

0.113 

u 

75 

0.200 

0.252 

0.284 

0.145 

.J 

77 

0.200 

0.282 

0.294 

0.162 

o 

>4 

22 

0.200 

0.295 

0.298 

0.170 

u 

85 

0.400 

0.335 

0.511 

0.193 

o 

87 

0.400 

0.32C 

0.506 

0.184 

z 

89 

0.400 

0.285 

0.495 

0.164 

15 

0.800 

0.415 

0.938 

0.239 

49 

0.800 

0.422 

0.940 

0.243 

70 

0.800 

0.380 

0.926 

0.219 

71 

0.800 

0.377 

0.925 

0.217 

108 

0.100 

0.187 

0.162 

0.107 

112 

0.100 

0.155 

0.151 

0.089 

115 

0.100 

0.180 

0.159 

0.103 

H 

114 

0.200 

0.253 

0.284 

0.146 

m 

132 

0.200 

0.300 

0.300 

0.173 

137 

0.200 

0.311 

0.303 

0.179 

H 

< 

133 

0.400 

0.412 

0.537 

0.237 

a 

136 

0.400 

0.425 

0.541 

0.245 

w 

137 

0.400 

0.425 

0.541 

0.245 

u 

130 

0.800 

0.390 

0.930 

0.225 

>4 

u 

134 

0.800 

0*440 

0.946 

0.254 

135 

0.800 

0.451 

0.950 

0.260 

138 

0.800 

0.250 

0.883 

0.144 

139 

0.800 

0*223 

0.875 

0.129 

107 

0.100 

0.170 

0.156 

0.098 

109 

0.100 

0.190 

0.163 

0.109 

113 

0.100 

0.165 

0.155 

0.095 

in 

105 

0.200 

0.280 

0.293 

0.161 

110 

0.200 

0.230 

0.276 

0.132 

H 

< 

123 

0.200 

0.270 

0.290 

0.155 

118 

0.400 

0.267 

0.489 

0.154 

Q 

w 

120 

0.400 

0.365 

0.521 

0.210 

•J 

121 

0.400 

0.363 

0.521 

0.209 

>4 

111 

0.800 

0.245 

0.881 

0.141 

u 

126 

0.800 

0.360 

0.919 

0.207 

128 

0.800 

0.410 

0.936 

0.236 

131 

0.800 

0.460 

0.953 

0.265 

TABLE  V 


INITIAL  SHEAR  MODULI 
(All  Stresses  in  KSI) 
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APPENDIX 


User's  Guides  and  FORTRAN  Listings  of  UNAX2  and  PROP. 

The  two  FORTRAN  programs,  UNAX2  and  PROP,  which  compute 
all  stress  and  strain  quantities  in  uniaxial  strain  and 
proportional  loading  tests,  respectively,  are  listed  in  this 
appendix.  A  description  of  the  input  and  output  for  each 
program  is  given.  All  the  computed  results  in  this  report 
were  obtained  by  execution  of  the  two  programs  on  the 
IBM  1130  computer  located  in  the  Weidlinger  office. 

UNAX2  utilizes  four  arithmetic  statement  functions, 

CFKLD(EE)  ,  CPKUN(P),  GLD(SQJ2,P)  and  GUN(SQJ?.,P)  which 

correspond  to  Eqs .  (3),  (4),  (5)  or  (9),  and  (6)  or  (10), 

respectively,  and  numerical  integration.  If  the  model 

description  were  to  change,  only  these  four  statements  need 

be  altered.  The  axial  stress  SIG  is  incremented  DSIG  and  a 

predictor  corrector  approach  is  used  to  find  the  strain  and 

remaining  stress  quantities  SIG3,  SI  and  P.  The  bulk  and 

shear  moduli,  K  and  G,  are  point  functions  at  the  new  stress 

and  strain,  as  are  the  (constrained)  tangent  and  secant 

4 

moduli,  EMTAN  and  EMSEC,  which  are  computed  from  K  +  —  G 
and  SIG/EPS.  The  local  wave  velocity  V  and  local  Poisson's 
ratio  NU  are  also  found  from  K  and  G.  The  typical  line  of 
output  consists  of  all  the  variables  mentioned  above.  In 
addition,  a  complete  heading  listing  all  input  parameters  and 
all  "useful"  constants  is  given.  The  program  loads,  unloads 
and  then  reloads  to  SIGL,  SIGU  and  SIGR,  respectively. 
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During  Initial  loading  the  Inflection  point  is  recorded 
and  the  information  Is  printed  at  the  end  of  the  loading 
phase.  To  reduce  the  amount  of  output  two  quantities  ISTOR 
and  IPRIH  are  introduced.  The  computer  stores  the  complete 
information  each  ISTORth  DSIC  on  the  d^sk  for  later  processing. 
Of  those  stored,  each  IPRIKth  is  printed  on  line. 


The  input  is  as  follows: 


Cards  1-9  FORMAT  (15A4) 


Card  10  FORMAT  (8E10.0) 


Card  11  FORMAT  (8E10.0) 


TITLE 

Alphanumeric  description  in 
columns  1-60  of  the  particular 
run  and  the  model .  The  first, 
printed  in  the  middle  of  the 
page,  usually  contains  identi¬ 
fication  of  the  particular  run. 
The  remaining  8,  printed  on 
alternate  sides  contain,  for 
example,  descriptions  of  K 
and  G. 

Values  of  the  material  para¬ 
meters  Ko  ,  KL  ,  K2  ,  Kqu  , 

K1U  *  V0  ’  ^1  *  Y1  (vo  is  the 
initial  Poisson's  ratio). 

Values  of  the  material  para¬ 
meters  y2  ,  Gqu  ,  ,  Ylw  » 

Pq  (Pq  is  the  initial  density) 


TITLE 


Card  12  FORMAT  (3A4,  E12.0)  Units  of  stress,  density  and 

velocity,  and  a  conversion 
factor  such  that 

CS  "  /G7P  *  *CONV  where  G.  p 
and  Cg  are  in  the  units  given 
e.g. ,  KSI ,  PCF ,  FPS  and 
2151.985. 

Card  13  FORMAT  (4E10.0,  215)  Values  of  program  controls 

SIGL ,  SIGU,  SIGR,  DSIC,  ISTOR 
IPRIN. 

Program  PROP  is  written  in  sections  for  loading,  un¬ 
loading,  reloading  of,  separately,  volumetric  and  deviatoric 
stresses.  Expressions  for  stress  and  strain,  Eqs.  (15),  (42), 
(43),  (48),  (49)  are  evaluated  explicitly  in  the  following 
manner . 

On  loading,  the  program  computes  stresses  at  discrete 
volumetric  strain  points.  The  strains  are  incremented  by  an 
amount  computed  in  the  program  as  the  volumetric  strain  at 
SlCLM(l)  divided  by  EENUM  -  the  input  number  of  strain  in¬ 
crements.  On  unloading  and  reloading  the  program  automatically 
shifts  to  stress  Increments  computed  as  the  difference  between 
the  pressure  at  the  loading  and  unloading  limits,  divided  by 
SGNUM  ■  the  input  number  of  stress  increments.  The  reason  for 
switching  from  strain  to  stress  increments  should  be  apparent 
from  Inspection  of  the  formulas  in  the  body  of  the  report 
appropriate  to  loading,  unloading  and  reloading. 
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In  Che  program  allowance  is  made  for  an  initial  set 
of  stresses  which  may  have  been  applied  to  the  test  sample. 
The  program  is  capable  of  handling  any  number  of  loading- 
unloading-reloading-loading  cycles  through  the  input  numbers 
SICLM(I)  which  are  values  of  specifying  points  at  which 
unloading,  reloading,  etc.,  are  to  occur. 

Data  is  input  to  PROP  according  to  the  following 


specifications . 

Card  #1  TITLE 

TITLE  (18) 

TITLE  (19) 

TITLE  (20) 

Card  #2  ANU 

(FORMAT  8E10.0)  RHOZ 

Q 

EENUM 

SGNUM 


Alphanumeric  title  in  columns 

1-68. 

Units  of  density,  i.e.,  PCF 
in  columns  69-72. 

Units  of  velocity,  i.e.,  FPS 
in  columns  73-76. 

Units  of  Btregg,  i.e.,  KSI 
in  columns  77-80. 

Poisson's  ratio  v. 

3 

Material  density  pQ  in  lb/ft  . 
Proportionality  constant 
(Q-0  gives  triaxial  test, 

Q-l  gives  Hydrostatic  test). 
Number  of  strain  increments 
to  be  taken  on  initial  loading. 
Number  of  stress  increments  to 
be  taken  on  unloading  and  re¬ 
loading  . 
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SIG1Z 

SIG2Z 

1  Initial  seating  stresses, 

I  VCT 

(J  j  ,  O  2  in  KS I  * 

Card  #3 

AKO 

Loading  constants 

(FORMAT 

8E10.0) 

AK1 

*  K^  ,  Ij  ,  ^2 

AK2 

bulk  modulus,  in  KSI. 

GMBR 

Loading  constants 

GAM1 

Yl  ,  .  Y2  for  shear 

GAM  2 

modulus,  in  KSI. 

Card  #4 

AKOU  ' 

Unloading  constants  K0(J  , 

(FORMAT 

8E10.0) 

AK1U 

for  bulk  modulus,  in  KSI. 

GMBRU  ] 

Unloading  constants  , 

GM1U  \ 

^  Gqu  for  shear  modulus. 

GOU 

in  KSI. 

Card  t 5 

NLIM 

Number  of  limit  loads 

(FORMAT 

15,  7E10.0) 

from  which  unloading 

1U 


'1U 


and  reloading  is  to  take  place 
SIGLM(I)  Values  of  stress  Oj  ,  in 

I»l,  NLIM  KSI ,  at  which  unloading  and 

reloading  occur. 


For  runs  with  multiple  sets  of  data,  insert  a  card  with  a 
in  column  1  between  the  data  sets 


U  V  KJ 


< 
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HYSTERMIC  MATERIALS 
UNIAXI.iL  STRAIN  TEST 
LOADING  -  UNLOADING  -  KELOAUING 
COMPILED  ON  JUNE  2S.  1968 
REAL  KZ.Kl.K2  .KZU.KlU  »NUZ .NU2U.M  31 ) «NU( 3i ) 

DIMENSION  SIG3(31).S1(31>.P(31) .tKKl H) .EMTAN(ii) .G<3i) .tMScC(Jl) • 
$  V( 3 1 ) t  TITLE (13*9) t SAVE (31*101 
EQUIVALENCE  ( SAVE ( 1 ♦ 1 1 • S 103 ( 1 1 > • i SAVE (1.* »  .Si ( 1 ) i . (SAVE ( 1 » 3 > *P( 1 ) > 

1  »<SAVE< 1*41 »FKK(1) »  .(SAVE* l.S) *K( l > ) « ( SAVE ( 1 .6 >  »G( 1 )  )  . (SAvtU *7) 

2  .EMTAN* 1) >• (SAVE* 1*8) »EMSEC( 1> ) • ( SAVE ( 1 # Y ) »V (i ) ) • ( SAVE l i » 1 J ) * 

3  NU(ll) 

C  DEFINE  FILES 

DEFINE  FILE  60 ( 15.30 «U« ITITL  > 

DEFINE  FILE  61(16. 60. U.ISIG3  > 

DEFINE  FILE  62( 16.60.U. IS1  > 

DEFINE  FILE  63(16*60. U. IP  ) 

DEFINE  FILE  64( 16.60 «U«  IEK.K  ) 

DEFINE  FILE  65<16.60iU.IK  ) 

DEFINE  FILE  66(16.60.U. IG  I 
DEFINE  FILE  67( 16*60*0. IMTAN  ) 

DEFINE  FILE  68( 16*60* U. IMSEC  > 

DEFINE  FILE  69( 16.60. U. IV  I 
DEFINE  FILE  70( 16.60. U. INU  » 

C  FUNCTIONS 

CPKLD(EE)  *  <2  ♦  EE**K1+K2*EE> 

CPKUN(P)  *  KZU  K1U»P 

AMINl(X.Y)  •  *5*( *1.-SIGN<1..X-Y) >*X+*1.+5IGN(1..X-Y) )*Y) 
GLD*SQJ2.P)=GZ+GAMH*SQJ2+AMINl(P.PC  MGAM1+GAM2*AM1M(p.PC> ) 

GUN ( S0J2  »  P ) *GZU+GAMBU*SQJ2+AM INI { P  »PC  >  * ( GAMiu+oAMZU*AMI N1  ( P  .PC ) I 
MOD  ( M .N )  «  M  -  M./N*N 

READ  (2.9991  TITLE 

C  READ  MATERIAL  PROPERTIES 

RE AD (2. 998)  KZ  *K1 *K2 .KZU.KlU.NUZ .GAME .GAM X .GAMZ .GZU .GAmBU .oAMIU 
$  .  RHOZ 

C  READ  UNITS  AND  CONVERSION  CONSTANT 

READ (2. 997)  PUNIT.DUNIT.VUNIT.VCONV 

C  READ  PROGRAM  CONTROLS 

READ (2.996)  SIGLtSIGU . S IGR .DSI G. I STOR » I  PR  IN 

C  COMPUTE  USEFUL  CONSTANTS 

HSQR3  «  SORT* 3.  )/  2. 

BE*A  «  2* * ( l.+NUZ ) / ( 1 »”?i*NUZ ) 

GZ  «  3**XZ/BETA 

PC  «  -GAM1/2./GAM2 

EZ  «  2»*( 1»+NUZ)*GZ 

EMZ  «  KZ  ♦  fc./3.*GZ 

CS  «  SORT (GZ/RHOZ)*VCONV 

CP  «  SORT ( FMZ/RHOZ ) *VCONV 

GAM2U*  GAM?*GAM1U/GAM1 

EMZU  «  KZU  ♦  4./3.+GZU 

CSU  *  SORT (  GZU/RHOZ )*VC0NV 
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CPU  «  SORT ( EMZU/RHOZ )  *VCONV 

BETAU  «  3  **KZU/GZU 

NUZU  «  0.5*( B£TAU“2. )/( BETAU  ♦  1.  > 

C  OUTPUT  TITLE.  INPUT.  AND  USEFUL  CONSTANTS  ON  DISK 

DO  25  NREC  «  1.9 

25  WRITE! 60  NREC  )  (  T I TLE ( J.NREC ) »J*1»15  ) 

WR I TE ! 60  10)  K2.K1.K2.KZU.K1U.NUZ.GAMB.GAM1.GAM2. 

S  GZU.GAMBU.GAM1U.RH0Z 

WRITF(6C  11)  PUN1 T .DUN I T . VUNI T .VCONV 
WRITE (60  1?)  S I GL • S I GU . S I GR  »  OS I G  » I  ST  OR  » I  PR  I N 

WRITE  1 60  13)  BETA.GZ.PC.EZ.EMZ »CS . CP. GAM2U.EMZU.CSU .CPU .BETAU. NUZU 

C  OUTPUT  INPUT 

WRITE! 3.995 )  TITLE 

WRITE (3.994)  KZ.K1.K2.GZ.PUN1T.GAMB.GAM1.GAM2.PUNIT.PC.EZ.EMZ. 

S  PUNIT.RH02.DUNIT.CP.CS.VUNIT.NUZ. BETA 
WRITE! 3.993)  KZU.PUNIT.KlU.GZU.PUNlT .GAMBU.GAM1U.GAM2U.PUNIT .EMZU. 
S  PUNIT.CPU.VUNIT.CSU.VUN  IT. BETAU. NUZU 
WRITE! 3.992)  SIGL.PUN 1 T.SIGU.PUNI T.S IGR .PUNI T .DSIG.PUNI T • 

S  ISTOR.  IPRIN 
WRITE! 3.991 )  PUNIT.VUNIT 
WRITE! 3.990  ) 

C  LOAD  TO  SIGL 

NMAX  «  SIGL/DSIG/FLOAT! ISTOR)  ♦  .001 
NRCMX  «  ( NMAX  +  29)/  30 
SIG  «  0. 

SIG3I1)  >  0. 

SKI)  >  0. 

P(l)  «  0. 

EKK(l)  «  0. 

K ( 1 )  «  KZ 
Gil)  >  GZ 
EMTAN(l)  ■  EMZ 
EMSEC(l)  «  EMZ 
VII)  «  CP 
NU(1)  =  NUZ 
EMIN  e  EMZ 
AKJ  «  K(l) 

GJ  *  G( 1) 

PJ  «  PIU 
S1J  «  Sill) 

EPSJ  =  EKK(l) 

EMJ  >  EMTAN(l) 

00  ISO  NREC«1» NRCMX 
DO  100  I  «1.10 
00  100  J  «  2.31 
100  SAVE ( J . I )  «  0.0 

NIST  «  30*(NREC-1)  ♦  1 
NLAST  ■  MINO(NMAX»30#NREC  ) 

00  150  NN  «  NIST  •  NLAST 
N  *  NN  -  NIST  +2 


00  130  J«1  * ISTOR 
DEPS  »  D5IG/EMJ 
OP  »  AKJ  /EMJ  *OSIG 
DS1  «  DS1G  -  OP 
EPSJ1  «  EPSJ  ♦  OEPS 
PJP1  »  PJ  ♦  DP 
S1JP1  «  SlJ  ♦  0S1 
AKJP1  «  CPKLDIEPSJ1/3.  ) 

GJP1  «  GLD(  HSQR3  *S1JP1*  PJP1  I 

EMJP1  «  AKJP1  +  GJP1/.75 

DEPS  «  0.5*(1./EMJ  ♦  1./EMJP1  )#DSIG 

DP  «  0.5*1 AKJ/EMJ  ♦  AKJP1/EMJP1  )*DSIG 

DS1  «  OSIG  -  DP 

SIG  «  SIG  +  DSIG 

PJ  «  PJ  DP 

S1J  «  S1J  ♦  DS1 

EPSJ  «  EPSJ  DEPS 

AKJ  «  CPKLD«EPSJ/3.  ) 

GJ  «  GLD(  HSQR3  *S1J«  PJ  ) 

EMJ  «  AKJ  ♦  GJ/.75 
IF  (  EMJ  -  EMIN)  110.130.130 
110  EMIN  «  EMJ 
SIGM  •  SIG 
PMIN  -  PJ 
EPSM  «  EPSJ 

VMIN  «  SORT ( EMIN/RHOZ )*VCONV 
130  CONTINUE 

SIG3IN)  >  PJ  -  S1J*0.5 
Sl(N)  *  S1J 
P ( N )  «  PJ 
EKK(N)  ■  EPSJ 
KIN)  *  AKJ 
GIN)  ■  GJ 
EMTANIN)  ■  EMJ 
EMSECIN)  «  SIG/EPSJ 

VIN)  «  SORT  I AMlNlI EMSEC IN) . EMJ ) /RHOZ )  *  VCONV 
NUIN)  ■  0.5*1 3»*K|N)  -2.*G(N) )/<3.*KIN)+GIN) > 

IF  I  MODINN.IPRIN) )  150.140.150 
140  WRITEI3.989  )  NN.SIG.  I  SAVEIN.J).  J>  1.10  ) 

150  CONTINUE 

JMAX  ■  MINO  I  N.  30  ) 

DO  160  NFILE  «  61.70 

160  WRITEINFIUE  NREC ) I  SAVE  I J.NFILE"60 > . J» 1 . JMAX ) 

DO  170  I  «1.10 

170  SAVE  1 1 . 1 )  *  SAVE!  H  *  I  ) 

180  CONTINUE 

WRITEI3.990  ) 

WRITE  I  3. 988)  SIGM .PUN  IT. PM  IN. PUN  IT. EPSM. EM  IN .PUN  IT  »VMl N  »VUN I T 
UNLOAD  TO  SIGU 

NMAXU  »  ISIGl  -  SIGU  ) /DSIG/FLOAT I  ISTOR)  +  .001 
NRCMU  *  I NMAXU  +  29)/30 
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E»SL  =  EPSJ 
PVAX  =  PJ 
K  {  1 J  *  CPK.UN(  PJ) 

G(l)  «  GUN(  HS0R3  *51  J»  PJ  » 

£VT AN ( 1 )  =  K(l)  ♦  G  < 1 ) /«  75 
EVSEC I 1 )  *  EvTAN(l) 

V < 1 1  «  SCRT(EMTAN(1)/RH0Z)*VC0NV 

NU  { 1 )  =  0.5*I3.*K'1)-Z.*G(1)I/(3.»KUI+0(1II 

AKJ  =  K ( 1 ) 

G  J  =  G  (  1 ) 

E"J  =  EMTAMD 
WR I TE ( 3  *9B7  1 
WR I Tf ( 3*990  1 
NREC  *  NRCMX 
NN  «  NMAX 

WRITE(3.9B9  )  NN.SIG.  <  SAVE(l.J).  J*  1.10  ) 
00  2B0  NR  »  l.NRCMU 
NREC  «  NREC  ♦  1 
DO  200  I  «i»ir 
DO  200  J  *2.3a 

200  SAVE ( J 1 1 )  *  0.0 

NLAST  =  M INO (NMAXU+1— ( NR-1 ) *30  »31 ) 

00  250  N*2. NLAST 
NN  «  NN-  1 
DO  230  J  =  l.ISTOR 
OEPS  -  OStG/EMJ 
OP  *  AKJ  /EMJ  *OSIG 
DS1  =  DSIG  -  OP 
EPSJ1  =  'IPSJ  -  OEPS 
PJP1  =  PJ  -  OP 
S1JP1  =  SlJ  -  OSl 
AKJP1  «  CPKUN ( P JP 1  ) 

IP(  S1JP1  )  210.215.215 

210  GJP1  «  GLO(-  HS0R3  *S1JP1»  PJP1  ) 

GO  TO  216 

215  GJP 1  *  GUN (  HSQR3  *S1JP1*  PJP1  ) 

216  GJP  1  =  -AMIMIO.O.  -GJP1  ) 

EWJP 1  =  AKJP1  ♦  GJP1/.75 

OEPS  *  0.5*<1./EVJ  1./EMJP1  )*DSIG 

DP  «  0.5MAKJ/EMJ  +  AKJP1/EMJP1  )*DSIG 

DSl  =  DSIG  -  DP 

SIG  =  SIG  -  DSIG 

PJ  «  PJ  -  DP 

SlJ  =  SlJ  -  DSl 

EPSJ  =  EPSJ  -  OEPS 

AKJ  =  CPKUN (  PJ  1 

IF  (  SlJ  )  220.220.225 

220  GJ  .  GLD<-  HSQR3  *S1J.  PJ  ) 

GO  TO  226 

225  GJ  «  GUN (  HS0R3  *S1J.  PJ  ) 

22 6  GJ  =  -AM INI (  0.0*  -GJ  ) 
tMJ  *  AKJ  +  GJ/.75 
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230  CONTINUE 

SIG3IN)  «  PJ  -  SI J*0.5 

Sl(N)  >  S1J 

P(NJ  «  PJ 

EKKINJ  «  EPSJ 

KIN)  >  AKJ 

GIN)  «  GJ 

EMTANIN)  ■  EMJ 

EMSECIN)  ■  ( SIGL  -SIGJ/IEPSL  -  EPSJ  ) 

VIN)  •  SQRTIEMJ/RHOZ)*VCONV 

NUIN)  •  0.5*1 3»*K(N)  -2.*GINi )/(3«*KCN)+GIN) ) 

IF  I  MODI NN* I  PR  IN) )  250.240.250 
240  WRITE  1 3*989  )  NN*SIG»  I  SAVEIN.,).  J«  1*10  ) 

250  CONTINUE 

JMAX  «  MINO  I  NLAST*  30  ) 

DO  260  NFILE  «  61»70 

260  WRITEINFILE  NREC) (SAVEI J»NFILE-60) .J-l.JMAX) 

DO  270  I  •  1*10 

270  SAVE  1 1 »  I )  >  SAVEI  NLAST* I  ) 

280  CONTINUE 

WRITE(3«990  ) 

C  RELOAD  TO  SIGR 

NMAXR  ■  I SIGR  -  SIGU  ) /DSIG/FLOAT I ISTOR )  ♦  .001 
NRCMR  «  I NMAXR  ♦  291/30 
EPSU  »  EPSJ 
IF  I  S1J  >  300.305.305 
300  GUI  ■  GUNI-  HSQR3  *S1J.  PJ  ) 

GO  TO  310 

305  Gil)  ■  GLDI  HS0R3  *S1J.  PJ  ) 

310  Gil)  ■  -AMIN1I  0.0.  -GIU  ) 

EMTANI1)  ■  Ml)  +  Gil)/. 75 

EMSEcm  ■  emtanid 

VI 1)  ■  SQRTIEMTANU  J/RHOZ)*VCONV 

Null)  •  0.5*13. *K<1)-2.*G<1) )/l3.*KI U+GIl) ) 

GJ  >  Gil) 

EMJ  ■  EMTANIl) 

WRITEI3.986  ) 

WRITEI3.990  ) 

WRITE (3.989  )  NN.SIG.  I  SAVE(l.J).  J>  1.10  ) 

DO  380  NR  ■  1. NRCMR 
NREC  •  NREC  +  1 
DO  312  I  •  1.10 
DO  312  J  ■  2.31 
312  SAVE(J.I)  >  0.0 

NLAST  ■  M INOINMAXR+i- (NR-1 ) *30 . 31 ) 

DO  350  N-2.NLAST 
NN  ■  NN+  1 
DO  330  J  ■  1. ISTOR 
DEPS  ■  DSIG/EMJ 
DP  ■  AKJ  /EMJ  *DSIG 
DS1  «  DSIG  -  DP 
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EPSJl  =  EPSJ  ♦  DEPS 

PJP1  «  PJ  ♦  DP 

S1JP1  «  SlJ  ♦  OS1 

IF  (  PJP)  -  PM AX  I  315*315*316 

315  AKJP1  =  CPXUNlPjPl  1 
GO  TO  ’17 

316  A< JP 1  «  CPXLDI EPSJl/3*  > 

317  IF  (  S1JP1  )  320  *320*321 

320  GJP1  =  GUM-  HSOR3  *S1JP1»  PJP1  ) 

GO  TO  322 

321  GJP 1  «  GLDI  HSQR3  *S1JP1*  PJP1  ) 

32?  GJP1  «  -AMINKO.O*  -GJP1  ) 

EMJP1  «  AX JP 1  ♦  GJP1/.75 

DEPS  *  0.5*1 l./EMJ  ♦  1./EMJP1  )*DSIG 

DP  «  0.5*<AKJ/EMJ  ♦  AX JP1/EMJP1  )*DSIG 

0S1  =  DSIG  -  DP 

SIG  *  SIG  ♦  DSIG 

PJ  ■  PJ  ♦  DP 

SlJ  *  SlJ  ♦  DS1 

EPSJ  *  EPSJ  ♦  DEPS 

IF  (  PJ  -  PMAX  )  323*324*324 

323  A< J  *  CPXUNl  PJ  ) 

GO  TO  325 

324  AX J  =  CPKLDl  EPSJ/  3.  ) 

325  IF  (  SlJ  )  327.32P.32P 

327  GJ  «  G'JNI-  HSQR3  *S1J*  PJ  ) 

GO  TO  329 

32P  GJ  «  GLD(  HSQR3  *S1J.  PJ  ) 

329  GJ  «  -AMINlt  0.0*  -GJ  > 

EMJ  >  AX J  ♦  GJ/ *75 

330  CONTINUE 

SIG3IN)  *  PJ  -  S1J*0.5 

S1(N)  =  SlJ 

PIN)  ■  PJ 

EXX(N)  ■  EPSJ 

X  ( N  )  «  AX  J 

GIN)  «  GJ 

EMTANIN)  ■  EMJ 

EMSECIN)  «  I S I GU  -S IG ) / I fcPSU  -  EPSJ  ) 

VIN)  «  SQRTIEMJ/RH0Z)*VC0NV 

NUIN)  «  0*5*1 3**XIN )  -?.*GIN) )  /  I  3  »*X  I  N  >*G  I  N  > ) 

IF  l  MOD(NN»IPRIN) )  350*340*350 
340  WR I TE 1 3  *989  I  NN*SIG*  I  SAVE(N*J)*  J»  1*10  ) 

350  CONTINUE 

JMAX  «  MINO  I  NLAST  *  3C  ) 

DO  370  NFILE  ■  61.  70 

WRITE  INFILE  NREC ) ( SAVE  I J *NF ILE-60 ) * J» 1  * JMAX ) 

370  SAVE! l.NF ILE-60)  *  SAVE  I JMAX* 1 *NF ILE-60  ) 

3«0  CONTINUE 

WRITE! 3.989  )  NN  *S IG . I  SAVE  I NLAST . J ) *  J-1.10  ) 
WRITE! 3*990  ) 

WR I TE l 60  14  )  NMAX«NRCMX.NMAXU*NRCMU»NMAXR*NRCMR* 
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S  SIGM.PMIN.EPSM.EMIN* VM1N 

C  FORMATS 

999  FORMATt  15A4  ) 

998  FORMAT (  8F10.2  ) 

997  FORMAT (  3A4,  E12.2  ) 

996  FORMAT(  4E10.2.  215  ) 

995  FORMATt  1H1  30X  15A4/  (IX  30A4  )  I 

994  FORMAT (31H0MATER I AL  PROPERTIES  -  LOADING  /  4H  KZ-F10.5.5H.  K.1* 

1  F10.3.  5H»  K2-F10.2.5H.  GZ*F10.5 .A4.7H.  GAMB*F10.4»7h »  6AM1* 

2  F10.4,  7H»  GAM2*  F10.5»2H  /  A4/  4h  PC-F8.3.  5H»  EZ»  F10.5. 

3  5H»  MZ-F10.5.A4.  7H.  RHOZ*  F7.3.A4.5H.  CP*  F8.3.  5H»  CS*  F8.3* 

4  A4»  6H»  NUZ*  F7.4.  7H.  BETA*  F7.3/  ) 

993  FORMATt  33H  MATERIAL  PROPERTIES  -  UNLOADING  /  5H  KZU*  F10.5»A4» 

1  6H»  XiU*  F10.4.  6H»  GZU*  F10.5.A4*  8H.  GAMBU*  F10.4.8H.  GAM1U- 

2  F10.4,  8H»  GAM2U-  F10.5.A4  /  5H  MZU*  F10.5»A4.  6H»  CPU*  F8.3i 

3  A4t  6H»  CSU*  F8.3»A4»  8H»  BETAU*  F8.3»7H»  NUZU*  F7.4/  ) 

992  FORMAT (  18H  PROGRAM  CONTROLS  /  6H  SIGL*F8.4.A4»  7H»  SIGU*  F8.4. 

1  A4t  7Hi  SIGR*  F8«4.A4*  7H»  DSIG*  F9.6.A4,  8H»  ISTOR*  I3» 

2  8H»  IPRIN*  13  /  ) 

991  FORMAT (  8H  LOADING  20X  20HALL  STRESSES  ARE  IN  A4» 

5  16H*  VELOCITIES  IN  A4  ) 

990  FORMATt /3H  N  6X  4HSIG1  6X  4HSIG3  7X  2HS1  9X  1HP  7X  4HEPS1  10X 
S  1HK  10X  1HG  9X  4HMTAN  6X  7HSIG/EPS  5X  1HV  7X  2HNU  /  ) 

989  FORMATt  14.  F10.3#  3FH.4.F10.5 »  4F11.3.F10.1  »F9.4  ) 

988  FORMATt  20H  POINT  OF  INFLECTION  /  7H  SIGM  *  F10.4.A4.5X  6HPMIN  * 

1  F10.4.A4.  5X  6HEPSM  •  F10.6*  5x  6HEMIN  •  F10.3.A4.SX  6HVMIN  * 

2  F10.1.  A4/  t 

987  FORMATt  10H  UNLOADING  ) 

986  FORMAT'  10H  RELOADING  I 
PAUSE 
CALL  EXIT 
END 
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C  PROPORTIONAL  LOAD  TEST 

C  SIGMA  2  DOT  *  0  *  SIGMA  1  DOT 

C  0*0  GIVES  TRIAXIAL  TEST 

C  0*1  GIVES  HYDROSTATIC  STRESS  TEST 

C 

DIMENSION  SIGLM<6) tTITLE(20)  .DP(6) 

C 

DATA  IASTE/  *  / 

I  N=2 
NOOT*3 

S3T3*S0RT ( 3 • ) 

S0T03*SQT3/3« 

C 

C  INPUT 

1000  CONTINUE 

READ! IN»910)  title 

10  READ! IN. ROD  ANU.RHCZ .0. EENUM.SGNUM.S I G1Z .SI G2Z 
C  LOADING  MODULI 

READt IN. 901)  AK0.AK1.AK2 .GM8R.GAM1 .GAM 2 
C  UNLOADING  MODULI 

READ < IN. 901 )  AK0U.AK1U.GMBRU.GM1U.G0U 
READ!  IN. 90?)  NL  I V  .  (  SI  GLMU)  .  1*1  .NL  IM  ) 

C 

C  NECESSARY  CONSTANTS 

LC YCL  *  1 
LV0L=1 
LDEV* 1 
LG0UT=1 
I 5AVE*0 
IPC-1 
EM0*0. 

0NTH*l./3. 

ONM IQ* 1 «-0 
0NP20*l.+2.*0 

PZERO* ( SI G1Z+2.*SIG2Z)*0NTH 
TOGM2*2.*GAM? 

SINIT*0*SIG1Z-SIG2Z 
CF* 1 • /0NP20 
CF 1 =ONM IQ*CF 
CF2*S0T03/CF1 
CF3*GMBR*CF2+GAMl 
CFA*GM1U/GAM1#GAM2 
G0=1.5*AK0*( l.-2.*ANU) / ( l.+ANU) 
GCPR*G0'.GMBR*S0T3*CF*SINIT 
G1BPL*SQT3*CF1*GMBR+GAM1 
G1PRU*SQT3*CF1*GMBRU+GM1U 
G1PCL*S0T3*CF1*GMBR 
G  1PCU*S0T3*CF  DGMBRU 
G?BRL*GAM2 
G2RRU=CF4 

RAD«SORT(G1BRL*G1BRL-4.*GOBR*G2BRL) 

I F  t  ABS  <  GAM? ) -# IE-3 )  11.11.12 


U  U  U  KJ  V  KJ 
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11  PLIM*-G0BR/G1RRL 
PC-l.E+6 

GO  TO  13 

12  PL IM«— ! G1BRL+RAD) /T0GM2 
PC—GAM1/T0GM2 
PC1».5*PC*GAM1 
IFIPlIM-PC)  13*13*1220 

1220  IFIO-1.)  1260*1250*1260 
1250  PLIM«1.E*6 
GO  TO  13 

1260  PLIM—  (G0BR+PCD/G1PCL 

13  CS«2153.3*SQRT ( G0/RH02 ) 

CP»CS*SORT(2.*< l.-ANU)/! l.-2.*ANU) ) 

SLIM1«3.*CF*IPLIM-PZFR0)+SIG1Z 
ELSCO«9.*GO*AKO/!GO+3.*AKO+Q*<  2.*G0-3.*AK0 ) ) 

00  14  I-2.NLIM.2 

DP( I)«<SIGlM( D-SIGLM! I-l> )/SGNUM*ONP2Q*ONTH 

14  DP ( 1+1 )«-DP ( I ) 

DEE« ( ONP20*ONTH*  !  S IGLN  !  1  » -S I G1Z 1 / < 3. *AKO ) ) /EENUM 
A*1.5*AK1 
B*3.0*AK0 
O-PZERO 
TRIAL*PZERO/B 

E0-FQB IC ! AK2. A. B.C. TRIAL. .00001.0) 

OUTPUT  OF  INPUT 

WRITE! NOUT .904 )  ( T I TLE ( I > . I« 1 . 17) 

WRITE(N0UT*914)  AK0.AK1.AK2. GO. TITLfcUO) .GMBR.GAMl.GAM2. TITLE! 20) » 

1  PC.E'.SCO  »RH0Z.TITLEI18).CP.CS.TITLEU9). 

2  ANU 

WRITE! NOUT .923)  AKOU.T I TLE < 20 ) .AKlU.GOU.T ITLE ! 20 ) .GMBRU.GM1U.CF4. 

1  T  ITlE I  20 ) 

WRITE! NOUT .913 )  «SIGLM! I) . T ITLE ( 20 ) . I* 1 .NLIM ) 

WRITE! NOUT .915 )  S IG1Z »T I TLE ( 20 ) . S IG2Z »T I TLE ! 20 ) .PLlM.T I TLE ! 20 ) . to 
WRITE!  NOUT  .916 )  SLIM 
WRITE! NOUT . 941 )  TITLEI20) 

WRITE!NOUT.950) 

IS»STEP  NUMBER 

LDEV* 1  FOR  DEVISTORIC  LOAD.  2  FOR  DEVIATORIC  UNLOAD 
IPC«1  FOR  P  LESS  THAN  PC.  2  FOR  P  GREATER  THAN  PC 

LOADING  CURVE 
G1BR«G1BRL 
G2RR«G2BRL 
G1PC-G1PCL 
PSAVE«PZERO 
PSG«PZERO 
GINIT-SINIT 
EESAV«E0 
E1SV«0. 

SG1SV«SIG1? 

S1«SIG1Z-PZER0 
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: ^1SV«C. 

C  VFASU«ED  LONGITUDINAL  STRAIN  IF^uM  SEATING)  COMPUTED 

c 

IS  tt*tESAV-DEc. 

Pl\C=0. 

Ccn*2.*G2^R*PSAVE+Gl«R 
RAD=SQRT(  Ci1HR*G1RR-4.*G0BR*G2RR) 
c 

03  ICO  1*1*1 000 
!S=  I  SAVE*- 1 

C  VOLUME TR I C  PORTION  -  LOADING 

f ESAV=LE 
fcp«FF»DEE 

1H  P*EL*(3.*AK0+EE*<l*5*AIU+AtC2*EE)  ) 

; F ( p— PL i v )  19*19*500 

19  A<=AKO+EE*(  AK  !♦  AtC 2*E F  > 

3IG1  =  CF*(  ,.*P*?.«SINIT) 

C  TEST  FOR  v A x I V jv  VALUE  OF  STRESS 

DIF*SIG1-SIGlM(LCVCL) 

IFtAHSlDIF)-l.E-4)  20*20.22 

20  LGOUT=? 

GO  TO  2A 

22  irtniFl  2A.2A.23 

2  3  P=OnP2C*0NTH* t S IGLVI LCYCL )~SIG1Z ) +PZERO 
O-P 
TRI AL=EE 

EF  =  FOHIC( ax  2, At H.C* TRIAL **00001*0) 

GO  U.  IB 

2a  S  !G;  *G*(  S  IC.l-S  I G)  2  )+S  IG2Z 
C  DEV  I  A  TuR I T  PORTION 

S1NF'.’*S  Kil-P 
IFIIF-1)  35.35.2AOO 

2^0  I F ( SINE W* ( SINEW-SI >  )  2A15. 35*2405 
?*.~5  I  F  t  LDEV-1  )  35.35.25 
2  a  1 5  ! f ( l^F V-2 )  25*35.35 

25  PSr,  =  EESAV*(  3.*AK0+tESAV*(  1 . 5*AH1+A'C2«EESAV )  ) 

E 1SV=E1 

29  GD  TO  (30.32) »LDEV 

C  DZVTATORIC  UNLOAD 

30  L OE V  =  2 

GTPR=GC+Gv«RU*S0T3*CF*GINIT 

r.IBwsMPPj 

G2PR*G?BRU 

G1PC=G1PCU 

(O  TC  3A 

C  DEV  I  A TOR  I C  LOAD  OR  RFLOAD 

12  LDE V= 1 

GOBR  =  GO*-GVRR*SOT3*CF*GINIT 

G1«R*G1RRL 

(>?HR*G2PRL 

34  RA-'sSORTI  G1RR*G1HR“A.*G0BR*G2BR  ) 

CF6*2.*G2P.R*PSG+G1RR 


r>  n  r\ 
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WR ITE ( NOUT #902  >  LDEV 

35  IF(P-PC)  36*3500*3500 
3500  GO  TO  (3505.40) *IPC 
3505  IPO? 

P«PC 

PSG«PC 

SIG1»CF*(3.*PC+2.*SINIT) 

SIG2*0*(SIG1-SIG1Z)+SIG2Z 

S1NEW-SIG1-P 

C=-PC 

TRIAL*EE 

EE«FOBIC( AK2.A.8.C* TRIAL* .00001*0) 

AK*AK0-*EE*(  AK1+AK2*EE ) 

G=G0BR+PC*(G1BP+GAM2»PC) 

CF5*2**G2BR*PC+G1BR 

E1SV«CF1/RA0*AL0G( I CF5-RADHM  CF6+RAD ) / ( (CF5+KAD)*(CF6-RAD) ) I+E1SV 
WRITEINOUT  *961) 

GO  TO  50 

36  IF(ARS(G2BR)-.lE-3)  39.37.37 

37  CF5=2.*G2RR*P+G1BR 
G*G0BR+P* ( G1RR+GAM2«P ) 

El«CFl/RAD*AlOG( < CF5-RAD)* ( CF6+RAD) / ( lCF5+RAD)«(CF6-RAD>  » I+E1SV 
GO  TO  50 

39  G*G0RR+G1BR*P 

El«CFl/GlRR*ALOG(G/(GORR«GiRR*PSG) )*E1SV 
GO  TO  50 

40  G=G0RR+PC1+G1PC*P 

E1*CF1/G1PC*AL0G(G/(G0BR+PC1+G1PC*PSG> 1+E1SV 
50  EPS1-E1+EE-EO 
SI«S1NEW 

ETAN»9.*G*AK/(ONP2Q*G-*3«*AH*ONMIO) 

65  ANU»0NMIQ*ETAN*«5/G-1. 

EPS2=\.5*(FE-E0)-.5*EPS1 
SOI  F*S  IC1-S  ir,2 
E0IF-EP31-EPS2 

WRITE (NOUT. 959) IS.SIG1 . S IG2 * 5DI F *P *EE *ED I F *EPS2 * EPS1 . AK .G.tTAN . ANU 
GO  TO  ( 100» 110) .LGOUT 
100  CONTINUE 

110  CONTINUE 
LGOUT  a 1 

120  IF(LCYCL-NLIM)  125*501.501 

UNLOADING  AND  RELOADING  CYCLES 

125  GO  TO  ( 130,140. 130) *LVOL 
130  WRITE (NOUT .957) 

LV0L»2 
A<K0»AK0U 
AKK1-AK1U 
PST«P 
GO  TO  160 
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U0  WRITE!  NOUT  ,  956 ) 

LV0L*3 

140  SGlSV=SIol 
SG2SV*SIG 2 
:P1SV*£PS1 
E 1 SV*E 1 
LCYCl*LCYCl*1 
I  SAVE* IS 
PSO*P 
P5AvE*P 
EE5AV*E£ 

CF6=2.*G2BR  *PSAVE*G1RR 

SIMT«0*SfclSV-Sr.2SV 

GIMT*SINIT 

P=PSAVE-PINC 

D.nl*DP(LCYCL) 

C 

no  300  1*1,1000 
IS* I  SAVE*  I 

C  VOLUMETRIC  portion  unload  or  reload 
Pxp+DPl 

SIG1*CF*(3.*P+2.*SINIT) 

IF(SIGI)  170,180,180 
1 70  S I G 1  *  0  , 

F1*-2.*0NTH*SINIT 
P I NC*P 1“P 
P  =  P1 

180  AK=AK0U*AK1U*P 

SIG2*Q*(SIG1-SG1SV)*SG2SV 
IF( AK1U-.1E-3)  185*190,190 
185  EE*CNTH*f P-PSAVE1/AK0U+EESAV 
GO  TC  195 

IPO  Et=ONTM*ALOG<AK/(AKOU+AKlU*PSAvE> ) 2AKlU*EtSAV 
195  CONTT'JE 

C  DEVIATORIC  PORTION 
200  S1NEW*SIG1-P 

I F ( SINEW* t S1NEW~S1 ) I  205,230,202 
202  IF(LDEV-l)  230,230*209 
205  IFILDEV-2)  209*230,230 
209  PSC'.*P-DP1 
E 1SV*F 1 

212  GO  TO  (215*220) ,LDEV 
C  LDFv=2  FOR  DEVIATOR IC  UNLOAD 
215  LDEV*2 

gorp*go+gmrru*sot3*cf*g In  it 

GlRPsGlflRU 

G2HR*G28RU 

G1PC=G1PCU 

GO  TO  225 

C  LCEV* 1  FOR  DEVIATOR  I C  LOAD  OR  RELOAD 

220  LDFV* 1 

G0RR*G0*GVRR*SQT3*CF*GINIT 


r>  n  n 
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G1BR-G1RRL 

G2BR*G2BRL 

G1PC-G1PCL 

225  CF6«2.*G2BR*PSG+G1BR 

RAD-SORT! G1BR*G1BR-4.*G0BR*G2BR) . 

WRITE! NOUT *902 )  LDEV 

230  IF!  P-PC  )  231*247.247 

231  GO  TO  <235*2321 . 1PC 

232  IPOl 
P=PC 

S1G1«CF*<3.*P02.*S1NIT> 

S IG2-Q*!  SIG1-SIG1ZI+SIG2Z 

S1NEW-SIG1-P 

C—PC 

TRIAL-EE 

EE-FQBIC! AK2.A.B.C. TRIAL. *00001.0) 

AK-AKC+EE*!  AKl-f  AK2*EE  ) 

G*G0BR+PC*!G1BR+GAM2*PC) 

CF5»2.*G2BR*PC+G1BR 

PSG-PC 

E1SV»CF1/RAD*AL0G<  <CF5-RAD)*!CF6*RAD> /< !CF5+RAD>*!CF6-RAD> ) 1+E1SV 
WRITE! NOUT .961 ) 

GO  TO  250 

235  IF(ABS!G2BR)-.lE-3)  245.240.240 
240  CF5*2.*G2BR*P+G1BR 

G-GOBR+P*! G1BR+G2BR*P  J 

E1=CF1/RAD*ALCG< !CF5»RA0)*(CF6+RA0)/< <CF5+RAD)*!CF6-RAD> ) I-.E1SV 
GO  TO  250 
245  G«G0BR+G1HR*P 

ElsCFl/GlBR*ALOG(G/(G0BR^GlBR*PSG))+ElSV 
GO  TO  250 

247  G»G0BR+PC1+G1PC*P 

EleCFl/GlPC*ALOG!G/(GOBR-*-PCl+GlPC*P5G>  I+E1SV 
250  EPS1-E1+EE-E0 
SleSlNEW 

FTAN«9.*G*AK/!ONP2Q*G+3.*AH*ONMIO> 

ANU»0NMI0*ETAN*.5/G-1. 

EP52«1.5*(EE-E0)-.5*FPS1 
EDI F«EPSl"FpS2 
SDIF-SIG1-SIG2 

WRITE (NOUT .959) I S.SIG1 »SIG2 »SDI F *P»EE »ED I F »tPS2 »EPS1 * AK .G.ETAN .ANU 
GO  TO  (260.260.270) .LVOL 
C  VOLUMETRIC  UNLOAD  TEST 

260  IF (A8S(SIG1-SIGLM(LCVCL) I-l.E-4)  125.125.300 
C  VOLUMETRIC  RELOAD  TEST 

270  IF(ABS(SlGl-SlGLM(LCYCLn-l.E-4)  125.125.280 
280  IF (ABS(P-PST)”l.E-4)  350.350.300 
300  CONTINUE 


VOLUMETRIC  RELOADING  ON  VIRGIN  CURVE 
350  WR I TE( NOUT .958 ) 


uuu 
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360  ! SAVE*  I S 
E1SV-E1 
FfcSAV=t£ 

PSAVE«P 
PSG-PST 
GO  TO  15 

600  WR I TE ( NOUT « 960 ) 

C01  REAO  !  2  #900 ) ICONT 

IF! ICONT-IASTE)  502 •2000*502 
50?  CALL  EXIT 

FORMATS 

POO  FORMAT! Al) 

P01  FORMAT  IflElO.O! 

902  FORMAT! I6*7E10.0) 

903  FORMAT { 3E20»7) 

90<*  FORMAT  {  1H1«//17A4//J 
°10  FORMAT  ( 20  A4 ) 

913  FORMAT ( /10H  SIGMA  1  • 

1  12HL0AD  LIMIT  *F10.5#A4.5X.14HUNLOAD  LIMIT  »F10«5*A4.5X 

1  14HREL0AD  LIMIT  «F10. 5# A4/ I50X, F 10.5 * A4 *29X#F10.5 tA4/ ) ) 

914  FORMAT ( 31H  MATERIAL  PROPERTIES  -  LOADING  /  4H  KZ*F10.5#5H#  Kl» 

1  Flo •  3 1  5h.  K2*F10#2  «5H»  GZ»F10.6tA4,7H.  GAMB*Fl0.4,7H#  GAMl* 

2  P10.4,  7h »  GAM 2*  F10.5#2H  /  A4/  4H  PC*F8.3.  5H.  EZ«  F10.5# 

3  5h*  7H •  RHOZ*  F7.3»A4.5h.  CP*  F8.3.  5H.  CS*  F8.3* 

4  A4.6H.  NUZ*  F7.4/) 

915  rJRvAT !  /1X» 14HSIGMA  1  ZERO  *F8. 5 .A4.4X . 14HSIGMA  2  ZERO  »Ed.5,A4. 

1  4X.9HP  LIMIT  =F8.5*A4.4X.20HMEAN  STRAIN  E  ZERO  *F8.5  ) 

916  FORMAT ( 55X» 15HSIGMA  1  LIMIT  *  FB.5.A4//I 

923  FORMAT  (  33M  MATERIAL  PROPERTIES  -  UNLOADING  /  5H  KZU*  F10.t>»A4. 

1  6H  *  K1U*  F10  #4 1  6H  *  GZU*  F10.5«A4.  BHt  GAMBU*  F10.4t8H*  GAM1U* 

2  F10.4.PH.  GAM2U*  F10.5*A4  /) 

041  FORMAT! 30X.20HALL  STRESSES  ARE  IN  A4/J 

950  FORMAT! /23X.4HSIG1*  HX*4MMEAN*5X»4HEPSl*4X*3HLAT»4Xt4HL0NG*25Xt4HT 
1ANG/25X.1H-.11X»6HSTRAIN#6X»1H-*1X*2!2X#6HSTRAIN) *24X»3HM0D/3X»1HN 
2.3Xt4HSIGl,4X.4HSIG3,4X*4HSIG3t6X*lHP»3X*6HEIOC/3..4X.4HEPS3*3Xt4HE 
3PS3 «4X #4MEPS1 *8X * 1HK  »8X • 1HG»  7X *4HMTAN  » 3X  »2HNU//8H  LOADING) 

956  FORMAT!  1CM  RELOADING  ) 

°57  FORMAT!  1CH  UNLOADING  ) 

958  FORMAT ! 2RH  RELOAD  ALONG  ORIGINAL  CURVE) 

959  FORMAT! 14 1 4F8 #4 t4F8 • 5 # 3F9 . 3  * F7 #4 ) 

960  F0RMAT14X.39HLIMIT  LOAD  REACHED#  NO  FURTHER  LOADING) 

961  FORMAT ( 5X  # 14HP  =  P  CRITICAL) 

END 
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The  present  report  extends  the  combined  variable  moduli  model 
Introduced  in  the  previous  report,. Ref.  (lj.  This  more  genera,  model 
is  defined,  conditions  are  set  on  the  various  parameters,  and  the  model 
behavior  in  uniaxial  strain,  triaxial  compression  and  proportional 
loading  tests  is  discussed. 

The  major  portion  of  the  present  report  deals  with  the  procedures 
used  to  fit  the  current  model,  including  the  loading  and  unloading,  to 
a  rather  complete  set  of  laboratory  data  for  McCormick  Ranch  Sand. 
Actually,  four  different  fits  are  described,  one  of  them,  Uniax-Triax  I 
in  some  detail.  The  theoretical  and  experimental  results  are  compared 
and  with  one  fit,  Uniax-Triax  II,  excellent  agreement  is  found  for 
uniaxial  strain,  triaxial  compression  and  proportional  loading  tests. 

Finally,  recommendations  are  made  concerning  reloading  in  shear. 
User's  guides  and  FORTRAN  listings  of  the  two  programs,  L*NAX2  and  PROP, 
used  to  compute  uniaxial  strain  and  proportional  loading  (including 
triaxial  compression;  tests  are  given  in  an  appendix. 
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